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Abstract 

Within the framework of QCD factorization (QCDF), we consider two different types of power correc- 
tion effects in order to resolve the CP puzzles and rate deficit problems with penguin-dominated two-body 
decays of B mesons and color-suppressed tree-dominated Ti^n^ and p^n^ modes: penguin annihilation 
and soft coiTcctions to the color-suppressed tree amplitude. We emphasize that the electroweak penguin 
solution to the B Kn CP puzzle via New Physics is irrelevant for solving the CP and rate puzzles re- 
lated to tree-dominated decays. While some channels e.g. , k"^ need penguin an- 
nihilation to induce the correct magnitudes and signs for their CP violation, some other decays such as 
B — > K^n^,n^ri,K^ri and B^ K*^ri,7i'^n^ require the presence of both power corrections to account 
for the measured CP asymmetries. In general, QCDF predictions for the branching fractions and direct 
CP asymmetries of S — > PP, VP, VV decays are in good agreement with experiment. The predictions of 
pQCD and soft-coUinear effective theory are included for comparison. 
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I. INTRODUCTION 



Although the underlying dynamics for the hadronic B decays is extremely complicated, it is greatly 
simplified in the heavy quark limit. In the nii, °° limit, hadronic matrix elements can be expressed in terms 
of certain nonperturbative input quantities such as light cone distribution amplitudes and transition form 
factors. Consequently, the decay amplitudes of charmless two-body decays of B mesons can be described 
in terms of decay constants and form factors. However, the leading-order 1 /mi, predictions encounter three 
major difficulties: (i) the predicted branching fractions for penguin-dominated B — > PP,VP,VV decays are 
systematically below the measurements |[l|] and the rates for color-suppressed tree-dominated decays B^ — > 
pOjj.0 small, (ii) direct CP-violating asymmetries for B K^K^ , B K*^K^, B K^p^, 

S — > TT+TT^ and Bg K^n^ disagree with experiment in signs, and (iii) the transverse polarization fraction 
in penguin-dominated charmless B VV decays is predicted to be very small, while experimentally it is 
comparable to the longitudinal polarization one. All these indicate the necessity of going beyond zeroth 
I /nib power expansion. 

In the QCD factorization (QCDF) approach power corrections often involve endpoint divergences. 
For example, the hard spectator scattering diagram at twist-3 order is power suppressed and posses soft 
and coUinear divergences arising from the soft spectator quark and the 1 / mi, annihilation amplitude has 
endpoint divergences even at twist-2 level. Since the treatment of endpoint divergences is model dependent, 
subleading power corrections generally can be studied only in a phenomenological way. Therefore, l/mt 
power suppressed effects are generally nonperturbative in nature and hence not calculable by the perturbative 
method. 

As a first step, let us consider power corrections to the QCD penguin amplitude of the B PP decay 
which has the generic expression 

P = Psd + Pld, 

= App[Xi,{a'l + r^al) + Xc{al + r^al)] + I /mh corrections, (1.1) 

where A^^' = VptV*^ with q = s,d, a4,6 ai^e the effective parameters to be defined below and is a chi- 
ral factor of order unity. Strictly speaking, the penguin contributions associated with the chiral factor 
are formerly l/m^, suppressed but chirally enhanced. Since they are of order l/m^ numerically, their ef- 
fects are included in the zeroth order calculation. Possible power corrections to penguin amplitudes include 
long-distance charming penguins, final-state interactions and penguin annihilation characterized by the pa- 
rameters P^''^. Because of possible "double counting" problems, one should not take into account all power 
correction effects simultaneously. As we shall see below in Sec. IV.B, CP violation of K^n^ and Tl^n^ 
arise from the interference between the tree amplitude Xlf-^ay and the penguin amplitude Xc^\a'^^ + ''^^6) 
with q = s for the former and q = dfor the latter. The short-distance contribution to a'^ + r^ag will yield a 
positive Acp(^^7r^) and a negative Acp(?r"'";r~). Both are wrong in signs when confronted with experiment. 
In the so-called "S4" scenario of QCDF |[l|l, power corrections to the penguin annihilation topology charac- 
terized by AujSj + Ac/33 are added to Eq. (11.11) . By adjusting the magnitude and phase of jSs in this scenario, 
all the above-mentioned discrepancies except for the rate deficit problem with the decays ff' — > 7C^7t^,p^7t^ 
can be resolved. 

However, a scrutiny of the QCDF predictions reveals more puzzles in the regard of direct CP violation. 
While the signs of CP asymmetries in K^7i^,K^p^ modes are flipped to the right ones in the presence 
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of power corrections from penguin annihilation, the signs of Acp in B —>■ K^n^, K rj, n^rj and 
71^71^, K*'^ri will also get reversed in such a way that they disagree with experiment. In other words, in 
the heavy quark limit the CP asymmetries of these five modes have the right signs when compai^ed with 
experiment. 

The so-called B — > Kn CP-puzzle is related to the difference of CP asymmetries of B K^n^ and 
5° K^K^. This can be illustrated by considering the decay amplitudes of S ^ Kn in terms of topological 
diagrams 

AiB'^^K-n^) = p' + T' + ^P^^ + P'^, 

AiB'^K'n') = ^(P'-C'-P^w-^^V + ^;), (1-2) 
AiB-^K^'n-) = P'-^Pi[^+A' + P[, 
A{B-^K-k'') = -L(p' + r' + C' + P^w + /EW+A' + /';), 

where T , C, E, A, Pew and P^^^ are color-allowed tree, color-suppressed tree, W-exchange, W-annihilation, 
color-allowed and color-suppressed electroweak penguin amplitudes, respectively, and Pa is the penguin- 
induced weak annihilation amphtude. We use unprimed and primed symbols to denote AS' = and jA^I = 1 
transitions, respectively. We notice that if C', P^^y and A' we. negligible compared with T' , it is cleai^ from 
Eq. (11.21 ) that the decay amplitudes of K^n^ and K^n^ will be the same apart from a trivial factor of 
1/a/2. Hence, one will expect that Acp{K^ n^) ^ Acp{K^K^), while they differ by 5.3a experimentally, 
^Kn =Acp{K-nP)-Acp{K-n+)=Q.U%±Qm% 

The aforementioned direct CP puzzles indicate that it is necessary to consider subleading power correc- 
tions other than penguin annihilation. For example, the large power corrections due to P' cannot explain 
the AAkti puzzle as they contribute equally to both B K^T^' and B" K^n^. The additional power 
connection should have little effect on the decay rates of penguin-dominated decays but will manifest in the 
measurement of direct CP asymmetries. Note that all the "problematic" modes receive a contribution from 
cC) = cC) + P^lf. Since A{B- R-n^) oc t' + c' + p' and A(S'' K-K+) oc t' + p' with t' = T' + P^^ and 
p' = P' — ^Pew + ^A' consider this puzzle resolved, provided that c' /t' is of order 1.3 ~ 1.4 with a 

large negative phase i\c' /t'\ ~ 0.9 in the standard short-distance effective Hamiltonian approach). There are 
several possibilities for a large complex c': either a large complex C' or a large complex electrowe ak p enguin 
P^w or a combination of them. Various scenarios for accommodating large C' . Inn orPg^ 

1 12|,ll3|] have been proposed. To get a large complex C', one can resort to spectator scattering or final-state 
interactions (see discussions in Sec.3.E). However, the general consensus for a large complex P^^ is that 
one needs New Physics beyond the Standard Model because it is well known that Pg^ is essentially real in 
the SM as it does not carry a nontrivial strong phase il4\ . In principle, one cannot discriminate between 
these two possibilities in penguin-dominated decays as it is always the combination c' = C +P^^ that enters 
into the decay amplitude except for the decays involving r] and/or r]' in the final state where both c' and Pg^ 
present in the ampUtudes lllSn . Nevertheless, these two scenarios will lead to very distinct predictions for 
tree-dominated decays where /few ^ C. (In penguin-dominated decays, Pg^ is comparable to C' due to the 
fact that Ac^' » Ai^^.) The decay rates of Ti^n''\p^n^ will be substantially enhanced for a large C but 
remain intact for a large Pew- Since Pew <C C in tree-dominated channels, CP puzzles with n^r\ and tt";!" 
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cannot be resolved with a large Pew- Therefore, it is most likely that the color-suppressed tree amplitude is 
large and complex. In other words, the B Kn CP puzzle can be resolved without invoking New Physics. 
In this work we shall consider the possibility of having a large color-suppressed tree amplitude with a 



sizable strong phase relative to the color-allowed tree amplitude jig] 



C= [A„a^]sD + [V^]LD + FSIs + --- . (1.3) 

As will be discussed below, the long-distance contribution to a2 can come from the twist-3 effects in spec- 
tator rescattering, while an example of final-state rescattering contribution to C will be illustrated below. 

Note that our phenomenological study of power corrections to penguin annihilation and to color- 
suppressed tree topology is in the same spirit of S4 and S2 scenarios, respectively, considered by Beneke and 
Neubert [;1]. In the "large S2 scenario, the ratio 0:2/^1 is enhanced basically by having a smaller Xb and 
a smaller strange quark mass. It turns out that the CP asymmetries of ,K*^ r] ,7i^ 

have correct signs in S4 but not so in S2, whereas the signs oi Acp{K^ K^),Acp{K^ r\)^Acp{K^K^) in S2 (or 
in the heavy quark limit) agree with experiment but not in S4. In a sense, our study is a combination of S4 
and S2. However, there is a crucial difference between our work and Jil], namely, our ^2 is not only large 
in the magnitude but also has a large strong phase. As we shall see, a large and complex a2 is needed to 
account for all the remaining CP puzzles. 

It should be remarked that the aforementioned B-CP puzzles with the K^Ti^, K rj, n^rj, K*^ri, n^Ti^ 
modes also occur in the approach of soft-collinear effective theory (SCET) linil where the penguin annihila- 
tion effect in QCDF is replaced by the long-distance charming penguins. Owing to a different treatment of 
endpoint divergence in penguin annihilation diagrams, some of the CP puzzles do not occur in the approach 



of pQCD III8II . For example, pQCD predicts the right sign for CP asymmetries of 5" — > ;r"7r" and ^ 7i^r\ 
as we shall see below. In this work, we shall show that soft power correction to the color-suppressed tree 
amplitude will bring the signs of Acp back to the right track. As a bonus, the rates of S° 71*' Tl", p^Ti^ can 
be accommodated. 

In the past decade, nearly 100 charmless decays of Bi,j mesons have been observed at B factories with 
a statistical significance of at least four standard deviations (for a review, see 111911 ). Before moving to the 
era of LHCb and Super B factories in the next few years, it is timing to have an overview on charmless 
hadronic B decays to see what we have leai^ned from the fruitful experimental results obained by BaBar and 
Belle. In this work, we will update QCDF calculations and compare with experiment and other theoretical 
predictions. 

This work is organized as follows. We outline the QCDF framework in Sec. 2 and specify various 
input parameters, such as form factors, LCDAs and the pai^ameters for power corrections in Sec. 3. Then 
Buj PP,VP,VV decays are analyzed in details in Sees. 4, 5 and 6, respectively. Conclusions are given in 
Sec. 7. 

II. B DECAYS IN QCD FACTORIZATION 

Within the framework of QCD factorization the effective Hamiltonian matrix elements are written 
in the form 

{M,M2\J^,n\B) = % I Xi>'\M,M2\^J^P+^^^'^P\B), (2.1) 

V ^ p=u,c 
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where the superscript h denotes the hehcity of the final-state meson. For PP and VP final states, h = 
0. J^'^ '' describes contributions from naive factorization, vertex corrections, penguin contractions and 
spectator scattering expressed in terms of the flavor operators af ''\ while :^^gg contains annihilation topology 
amplitudes characterized by the annihilation operators bf '^\ Specifically 

^J' = aP{MiM2)5pu{ub)^_^(g) {qu)^_^+a^{MiM2)5pu{qb)^_^(g) {au)y_^ 

+ af (MiM2)£(#),_^ {q'q'),_,+al{MxM2)Y,{q'b\_^ ® {qq\_^ 

+ 4(MiM2) {q'q%^,+a'^{MiM2)Y,{-2){q'b),_, {qq'),^, (2.2) 

+ fl^(MiM2) ^eg{q'qX^,+aP{MiM2)Y,i-2)iqb),_, (E)^iqq'),^, 

3 3 
+ af(MiM2) ® -e^{q'q\_^+aPQ{MiM2)Y,{q'b),_, ® -eci{qq')y_„ 

where {qiq2)v±A =^iYi^{^^ Ys)^! and {q\q2)s±p = =t 75)^2 and the summation is over = u,d,s. The 
symbol indicates that the matrix elements of the operators in are to be evaluated in the factorized 
form. For the decays B — > PP,V P,VV , the relevant factorizable matrix elements are 

{P2\J^\0) (Pi \J'^ \B) = ifp, {ml - ml {m\), 

(V|7^|0)(P|/;|B) = IfvmBPcFfiml), 



{P\J^\0){V\J'\B) = IfpmBPcAf'imi) 



X 



{BVi,V2. 
h 



{V2\J^\0){Vy\j'^\B) = -ifv,m2 



(ei*-e2*)K + mvjAf'(m; 



i£i ■ Pb) {£2 ■ Pb 



[m 



V2> 



(ms + my. 



+ 't^vaj3t2 ^1 PbP\ 



(2.3) 



(niB + mvi] 

where we have followed the conventional definition for form factors 12011 . For B — > VP,PV amplitudes, we 
have applied the replacement mv{s* • ps) — > msPc with pc being the cm. momentum. The longitudinal 
{h = 0) and transverse (h = it) components of xjf^^'^^'^ are given by 



y(BVi,V2) 
^0 



ifv2 



2 „2 



X. 



{BVi,V2) 



2mvi 
-ifv^mBmv2 



{ml - < - ml){ms + my, )Af ' {q^) - -±!1mE^A 



l + ^!^Vf'(^^)T^^ 
ms J mB + mvj 



mg + wivi 
V^^'{q^) 



BV, 
2 



(2.4) 



The flavor operators af '' are basically the Wilson coefficients in conjunction with short-distance non- 
factorizable corrections such as vertex corrections and hai^d spectator interactions. In general, they have the 
expressions 



{M1M2 



N'l{M2) + - 



Vl^{M2) + —H'}{M,M2] 



+ P';^P{M2), (2.5) 



Nc J Nc 4k 

where / = I,-- - ,10, the upper (lower) signs apply when / is odd (even), c, are the Wilson coefficients, 
Cp = {N^ — 1)/ {2Nc) with Nc = 3, M2 is the emitted meson and Mi shares the same spectator quark with the 
B meson. The quantities V/'(M2) account for vertex corrections, H'1{M\M2) for hard spectator interactions 
with a hard gluon exchange between the emitted meson and the spectator quark of the B meson and Pi{M2) 
for penguin contractions. The expression of the quantities N'1{M2) reads 



N^{M2 




(2.6) 
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The weak annihilation contributions to the decay B —>■ M1M2 can be described in terms of the building 



blocks bf''^ and ^f^w 
Gf 



V p=u,c V ^ p=u,c i 



(2.7) 



The building blocks have the expressions yj 

^1 = ^^1^1' 



^3,EW 



J4,EW 



N} 



C3A'i+C5(A^+A{)+A^,C6A( 



C4A\ + C6A{ 



(2.8) 



Here for simplicity we have omitted the superscripts p and h in above expressions. The subscripts 1,2,3 of 
A'n-^ denote the annihilation amplitudes induced from (V -A){V -A), {V -A){V+A) and {S-P){S + P) 
operators, respectively, and the superscripts / and / refer to gluon emission from the initial and final-state 
quarks, respectively. Following yj] we choose the convention that Mi contains an antiquark from the weak 
vertex and M2 contains a quark from the weak vertex. 

For the explicit expressions of vertex, hard spectator corrections and annihilation contributions, the 
reader is referred to III, Q, llll] for details. The decay amplitudes of 5 ^ PP, VP are given in Appendix 
A of III] and can be easily generalized to B ^ VV (see 1 22] for explicit expressions of B — > VV amphtudes). 
In practice, it is more convenient to express the decay amplitudes in terms of the flavor operators cu^'' and 
the annihilation operators which are related to the coefficients a^'' and bf by 



a';{MiM2) 
(4{MiM2) 

a!^'P{MiM2) 



al'P{M\M2) 



a\{M,M2), 
4{M,M2), 
{ a/{MiM2 

\ 4''\MiM2)+a"/{MiM2) 
{ c^l''{M,M2) + rfaY{MyM2) 



-a,'\MxM2) 



for M\M2 
for M1M2 

for M1M2 

(M1M2) - (M1M2 ) for M1M2 



PP, VP, 
VV, PV, 



Ml h,P ( 



■ PP, PV, 
VP,VV, 



(2.9) 



a 



h,p 
3,EW 



(M1M2 



h,p 



{M1M2 



■ a 



'l'P{M,M2)+c^'/{M,M2) 



{M1M2 



for M\M2 
for M1M2 



PP, VP, 
VV, PV, 



a: 



h,p 



4^(MiM2) + rf a7(MiM2) 



M2 h,P I 



'4,EW 



(M1M2 



a';o"(MiM2)-rf a7(MiM2) 



Ml h,p I 



for M\M2 
for MyM2 



■ PP, PV, 
VP,VV, 



and 



ifBfMifM2 u> 
X(BMi,M2) ' ■ 



(2.10) 



j8f(MiM2) 

The order of the arguments of af {M1M2) and jSf (M1M2) is consistent with the order of the arguments of 
^{BMiMi) =^^^^^_ jj^e chiral factor is given by 

2mv f^iix) 



2mp 



m/,(jU)(m2+"n)(M)' 



mhiP-) fv 



(2.11) 
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The Wilson coefficients c,(/i) at various scales, = 4.4 GeV, 2.1 GeV, 1.45 GeV and 1 GeV are taken 
from 12311 . For the renormalization scale of the decay amplitude, we choose jj. = mb{mh). ^ However, as 
stressed in the hard spectator and annihilation contributions should be evaluated at the hard-coUineai^ 
scale Hh = a/mA/, with A/, « 500 MeV. 



m. INPUT PARAMETERS 
A. Form factors 

There exist many model calculations of form factors for B — > P, V transitions. For B — > P transitions, 
recent light-cone sum rule results for form factors at = are collected in Table U A small F^^ of order 
0.25 is also preferred by the measurement of B — > n^n^. It is more convenient to express the form factors 
for B — > T] transitions in terms of the flavor states qq = {uu + dd) /V2, ss and cc labeled by the T]^, T]., and 
T]^, respectively. Neglecting the small mixing with T]'', we have 

^Br, ^^B77,cos0, F^''' = F^''''sin0, (3.1) 

where Q is the tj^ — tj^ mixing angle defined by 

|t]) = cos0|t]^) -sine|T],), 

\^') = sin0|T]^) + cos0|T],), (3.2) 



with Q = (39.3 ± 1.0)° in the Feldmann-KroU-Stech mixing scheme ll25ll . From the sum rule results shown 
in Table U we obtain ^^^''"(0) = 0.296. The flavor -singlet contribution to the B — > T]'^ form factors is char- 
acterized by the parameter B|, a gluonic Gegenbauer moment. It appears that the singlet contribution to the 



form factor is small unless Bj assumes extreme values ~ 40 II29I1 . 

The B ^ n,K, rjq transition form factors to be used in this work are dispalyed in Table HIl We shall use 
the form factors determined from QCD sum rules for B ^ V transitions 1301]. 



' In principle, physics should be independent of the choice of /i, but in practice there exists some residual ji depen- 
dence in the truncated calculations. We have checked explicitly that the decay rates without annihilation are indeed 
essentially stable against ji. However, when penguin annihilation is turned on, it is sensitive to the choice of the 
renormalization scale because the penguin annihilation contribution characterized by the parameter is dominantly 



proportional to as{Hh)c(,i^ih) the hard-colHnear scale = ^juA/,. In our study ofB—>VV decays 12411 . we found 
that if the renormalization scale is chosen to be /i = mi,{mi,) /2 = 2.1 GeV, we cannot fit the branching ratios and 
polarization fractions simultaneously for both B ^ K*(^ and B K*p decays. In order to ensure the validity of the 
penguin-annihilation mechanism for describing B VV decays, we will confine ourselves to the renormaUzation 
scale jj. = mt,{mi,) in the ensuing study. 
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TABLE I: Form factors for B ^ P transitions obtained in the QCD sum rules with B| being the gluonic 
Gegenbauer moment. 



F,f(0) 


0.258 ±0.031 \26} 


26+° °'^ 


[27] 


f;r(0) 


0.331 ±0.041 \26} 


36+^-^^ 

"•-'°-0.04 


[2a 




0.229 ±0.024 ±0.011 [2^] 




0. 188 ± 0.002Bf ± 0.019 ± 0.009 [2^] 



B. Decay constants 



Decay constants of various vector mesons defined by 

{V{p,e)\q2Y^qi\0) = -ifvmve*, 



(3.3) 



are listed in Table |TI1 They are taken from 13 ih . For pseudoscalar mesons, we use /;t = 132 MeV and 



fx = 160 MeV. Decay constants /^(/j, /^(/) and defined by 



are also needed in calculations. For the decay constants and /^(/, , we shall use the values 
4 = 107MeV, /^ = -112MeV, /^^, = 89MeV, = 137 MeV 



(3.5) 



obtained in i25ll . As for /'^(/), a straightforward perturbative calculation gives 13211 



12m2 y/2' 



(3.6) 



C. LCDAs 



We next specify the light-cone distribution amplitudes (LCDAs) for pseudoscalar and vector mesons. 
The general expressions of twist-2 LCDAs are 



^p{x,ix) = 6x(l -x) 



^\{x,pl) = 6x(l-x) 



^\{x,ix) = 6x{l-x) 



i + j^j;,{pi)cl'^{2x-\) 

l + Y^al{^)cT{2x-l) 

OO 

l + Y^aj^-y{ll)cl'\2x-\ 



n=\ 



(3.7) 
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TABLE II: Input pai^ameters. The values of the scale dependent quantities fy{pi) and {^i-) are given for 
= 1 GeV. The values of Gegenbauer moments are taken from 13311 and Wolfenstein parameters from 1341] . 



Light vector mesons 


V 


/v(MeV) 


/i^(MeV) 










P 
CO 

</> 

K* 


216±3 
187 ±5 
215±5 
220 ±5 


165 ±9 
151±9 
186±9 
185 ±10 






0.03 ±0.02 


0.15 ±0.07 
0.15 ±0.07 
0.18 ±0.08 
0.11 ±0.09 






0.04 ±0.03 


0.14 ±0.06 
0.14 ±0.06 
0.14 ±0.07 
0.10 ±0.08 


Light pseudoscalar mesons 


flf 




af 







0.25 ±0.15 


0.06 ±0.03 


0.25 ±0.15 


B mesons 


B 


mB(GeV) 


Tb(ps) 


/B(MeV) 


AB(MeV) 


Bu 


5.279 


1.638 


210 ±20 


300 ±100 


Bd 


5.279 


1.525 


210 ±20 


300 ±100 


Bs 


5.366 


1.472 


230 ±20 


300 ±100 


Form factors at = 


Fo«^(0) 




Af^*(0) 


Af^*(0) 




0.35 ±0.04 


0.374 ±0.033 


0.292 ±0.028 


0.259 ±0.027 


0.411 ±0.033 






Af (0) 


Af(0) 


C(o) 


0.25 ±0.03 


0.303 ±0.029 


0.242 ±0.023 


0.221 ±0.023 


0.323 ±0.030 


C(o) 


A^(0) 


Af«'(0) 


Af'»(0) 


V,f-(0) 


0.296 ±0.028 


0.281 ±0.030 


0.2 19 ±0.024 


0.198 ±0.023 


0.293 ±0.029 


Quaik masses 


mi,{mh) /GeY 


mc{mh) /GeV 


pole / pole 

mi /ml 


m,(2.1 GeV) /GeV 


4.2 


0.91 


0.3 


0.095 ±0.020 


Wolfenstein parameters 


A 


X 


p 
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0.8116 


0.2252 


0.139 


0.341 


(67.81^1)° 



and twist-3 ones 



%{x) = 1, 



<I>fj(x) = 6x(l -x), 

00 

2x-l + £4'^(^)P„+i(2x-l) 



n=l 



(3.8) 
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where C„(x) and Pn{x) are the Gegenbauer and Legendre polynomials, respectively. When three-particle 
amplitudes are neglected, the twist-3 <I>v(x) can be expressed in terms of ^± 



du 



du. 



The normalization of LCDAs is 



[ dxf£>vix) = l, [ dx(£>y{x)=0. 
Jo JO 



(3.9) 



(3.10) 



Note that the Gegenbauer moments a] 
ing a strange quark. 

The integral of the B meson wave function is parameterized as |Ql 

"1 dp 



displayed in Table Ull taken from i33h are for the mesons contain- 



1 



-a>f(p) = — , 



(3.11) 



where 1 — p is the momentum fraction carried by the light spectator quark in the B meson. The study of 
hadronic B decays favors a smaller first inverse moment Xb- a value of 350 ± 150 MeV was employed in iQl 
and 200;^^° MeV in UM, though QCD sum rule and other studies prefer a larger ~ 460 MeV 1351]. We 
shall use Xb = 300 ± 100 MeV. 



For the running quark masses we shall use 1361 .1371] 

mfo(mfo) =4.2GeV, m/,(2.1GeV) = 4.94GeV, 

mc{mt)= 0.91 GeY, mc(2.1GeV) = 1.06GeV, 

m,(2.1GeV) =95MeV, m,(lGeV) = llSMeV, 
mrf(2.1GeV) = 5.0MeV, m„(2.1GeV) = 2.2MeV. 



mfo(lGeV) =6.34GeV, 
mc(lGeV) = 1.32GeV, 



(3.12) 



Note that the charm quark masses here are smaller than the one mc{mh) = 1.3 ± 0.2 GeV adopted in HI, 122 
and consistent with the high precision mass determination from lattice QCD 1381] : mc(3GeV) = 0.986 ± 
0.010 GeV and mc{mc) = 1.267 it 0.009 GeV (see also 1391] '). Among the quarks, the strange quark gives 
the major theoretical uncertainty to the decay amplitude. Hence, we will only consider the uncertainty in the 
strange quark mass given by ms{2.\GdV) = 95 it 20 MeV. Notice that for the one-loop penguin contribution, 
the relevant quaik mass is the pole mass rather than the current one 1401] . Since the penguin loop correction 



is governed by the ratio of the pole masses squared Si = [n-^^'^ /ni^^^'^Y since the pole mass is meaningful 
only for heavy quarks, we only need to consider the ratio of c and b quark pole masses given by Sc ~ (0.3)^. 



D. Penguin annihilation 

In the QCDF approach, the hadronic B decay amplitude receives contributions from tree, penguin, elec- 
troweak penguin and weak annihilation topologies. In the absence of 1 / power corrections except for 
the chiral enhanced penguin contributions, the leading QCDF predictions encounter three major difficulties 
as discussed in the Introduction. This implies the necessity of introducing l/m/, power corrections. Soft 
corrections due to penguin annihilation have been proposed to resolve the rate deficit problem for penguin- 
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dominated decays and the CP puzzle for —>■ K^Ti^. ^ However, the penguin annihilation amplitude 
involve troublesome endpoint divergences. Hence, subleading power corrections generally can be studied 
only in a phenomenological way. We shall follow [2] to model the endpoint divergence X = Jq dx/x in the 
annihilation and hard spectator scattering diagrams as 



XA=ln 



(3.13) 



with A/, being a typical scale of order 500 MeV, and Pa, (I>a being the unknown real parameters. 

A fit to the data of Bu,d PP,VP,PV and VV decays yields the values of Pa and 0a shown in Table [ 
Basically, it is very similar to the so-called "S4 scenario" presented in l^. Within the framework of QCDF, 
one cannot account for all charmless two-body B decay data by a universal set of Pa and 0a parameters. 
Since the penguin annihilation effects are different for B VP and B — > PV decays, 

t2' 



for M1M2 = VP and 



67ra.v 



■ 67rav 



-3rl 



3(xr-A+'L\+//({xr?_2X^ 



3rl{lxY-l){l 



2XY + ^ 



X' X 

2 



^/A(xYf-ixY^ 



X 



VP\ 



/MxYf- 



vVP 



(3.14) 



xr-A+'^)+ry,{{xrf-2x, 



PV 

A 



1/ 



6Kas 



3r^ ( {Xrf 



2Xr + 4- 



71" 



71" 



3r'^{2XY-m-Xr) + rl{2 {X^ f-X^ 

rVP 



rPV 



(3.15) 



PV, the parameters Z( and X^ 



are not necessarily the same. Indeed, a fit to the B —>■ VP, PV 
70° and p^ ~ 0.87, ~ "30° (see Table |III1). For the estimate of 



for M1M2 

decays yields pj^ 1.07, 

theoretical uncertainties, we shall assign an error of ±0.1 to Pa and ±20° to (pA- Note that penguin annihi- 
lation contributions to K(j) {K*(j)) are smaller than other PV (VV) modes. In general, penguin annihilation is 
dominated by or jSs through {S — P){S + P) interactions. 



TABLE III: The parameters Pa and <pA for penguin annihilation. The fitted Pa and (pA for B 
aie taken from 1241] . 



VV decays 



Mode Pa (j)A 


Mode Pa 0a 


B^PP 1.10-50° 
B^PV 0.87 -30° 
B^K*p 0.78 -43° 


B^VP 1.07 -70° 
B^Kcj) 0.70 -40° 
B^K*<p 0.65 -53° 



Besides the mechanisms of penguin annihilation, charming penguins and final-state rescattering, another possibility 
of solving the rate and CP puzzle for K ;r+ was advocated recently in 14 If by adding to the B Kn QCDF 
amplitude a real and an absorptive part with a strength 10% and 30% of the penguin amplitude, respectively. 
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E. Power corrections to ^2 



As pointed out in jlql . while the discrepancies between theory and experiment for the rates of penguin- 
dominated two-body decays of B mesons and direct CP asymmetries of B^i — > K^n^, B K^p'^ and 
Bd K'^K~ are resolved by the power corrections due to penguin annihilation, the signs of direct CP- 
violating effects in B —>■ K^7i^,B^ —>■ K rj and B° —>■ n^n^^ are flipped to the wrong ones when confronted 
with experiment. These new B-CP puzzles in QCDF can be explained by the subleading power corrections 
to the color-suppressed tree amplitudes due to hard spectator interactions and/or final-state interactions 
that yield not only correct signs for aforementioned CP asymmetries but also accommodate the observed 
Bd n^nP and p^K^ rates simultaneously. 



Following jlql . power corrections to the color-suppressed topology ai^e parametrized as 



ai -^a2{l+Pce" 



with the unknown parameters pc and 0c to be inferred from experiment. We shall use jlql 

PC -1.3, 0.8, 0, (/)c~-70°,-80°, 0, 



(3.16) 



(3.17) 



for B —>■ PP,VP,VV decays, respectively. This pattern that soft power corrections to a2 are large for PP 
modes, moderate for VP ones and very small for VV cases is consistent with the observation made in |Q] 
that soft power correction dominance is much larger for PP than VP and VV final states. It has been argued 
that this has to do with the special nature of the pion which is a qq bound state on the one hand and a nearly 
massless Nambu-Goldstone boson on the other hand [9]. 

What is the origin of power corrections to a2 ? There are two possible sources: hard spectator interactions 
and final-state interactions. From Eq. (13.181 ) we have the expression 

ci CpOts 



£?2 M1M2) = C2 + — + — V2(M2 + //2 M1M2 



+ £?2(MiM2)lD, 



for ^2- The hard spectator term 7/2(^1 ^2) reads 

H2{M,M2) = ? C dxdy ( ^^^^^)^^^^y) + ^'nM)^MM 



where x^^'^^' '^^) the factorizable amphtude for B M1M2, x ■■ 
amplitude 0„, ai^e divergent and can be parameterized as 



(3.18) 



(3.19) 



I —X. Power corrections from the twist-3 



X. 



H 



^^^=ln^^(l+p..'-^- 
y Ah 



(3.20) 



Since ci ~ and cg ~ (^(—1.3) in units of it is clear that hard spectator contributions to tz, are 
usually very small except for a2 and a\o. Indeed, there is a huge cancelation between the vertex and naive 
factorizable terms so that the real part of 02 is g overned by spectator interactions, while its imaginary part 
comes mainly from the vertex corrections i42i 1. The value of a2{Kn) w 0.51e^'^^° needed to solve the 
B Kn CP puzzle [see Eq. (14.41 )1 corresponds to p^ ~ 4.9 and 0// k, —IT . Therefore, there is no reason 
to restrict pu to the range < p// < 1. A sizable color-suppressed tree amphtude also can be induced via 
color-allowed decay B K rj' followed by the rescattering of K rj' into K^n^ as depicted in Fig. 1. 
Recall that among the 2-body B decays, B — > Kr]' has the largest branching fraction, of order 70 x 10^^. 
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u 




II. u 



FIG. 1: Contribution to the color-suppressed tree amplitude of B K^K from the weak decay B 
K rj' followed by the final-state rescattering of K rj' into K^n^. This has the same topology as the color- 
suppressed tree diagram. 



This final-state rescattering has the same topology as the color-suppressed tree diagram i43h . One of us 
(CKC) has studied the FSI effects through residual rescattering among PP states and resolved the B-CP 
puzzles 101. As stressed by Neubert sometime ago, in the presence of soft final-state interactions, there is no 
color suppression of C with respect to T 



441. 



Since the chiral factor for the vector meson is substantially smaller than for the pseudoscalar meson 
(typically, = ^(0.8) and = ^(0.2) at the hard-coUinear scale jLt/, = -y/Aftmi), one may argue that Eq. 
( 13.191 ) provides a natural explanation as to why the power connections to 02 is smaller when Mi is a vector 
meson, provided that soft corrections arise from spectator rescattering. Unfortunately, this is not the case. 
Numerically, we found that, for example, H{K*n) is comparable to H{K7i). This is due to the fact that 
jQdxr^(i>m{x)/{l -x) is equal toXnr^ for M = 7^ and approximated to 3{XH-2)r^ for M = V. 

We use NLO results for a2 in Eq. (13.161 ) as a benchmark to define the parameters pc and (pc- The NNLO 
calculations of spectator-scattering tree amplitudes and vertex corrections at order aj have been carried 
out in 14511 and 14611 . respectively. As pointed out in 1471. 14811 . a smaller value of Xb can enhance the hard 
spectator interaction and hence az substantially. For example, aiinn) ~ 0.375 — 0.076/ for = 200 MeV 
was found in However, the recent BaBar data on S ^ ytv \49] seems to imply a larger Ag (> 300 MeV 
at the 90% CL). While NNLO corrections can in principle push the magnitude of a2{7in) up to the order 
of 0.40 by lowe ring the value of the B meson parameter Ag, the strong phase of a2 relative to ai cannot be 
larger than 15° 14711 . In this work we reply on pc and (j)c to get a large magnitude and strong phase for 02- 



rV. B^PP DECAYS 

Effects of power corrections on penguin annihilation and the color-suppressed tree amplitude for some 
selective B PP decays are shown in Table UV] The implications will be discussed below. Branching 
fractions and CP asymmetries for all B — > PP decays are shown in Tables |V] and IVIII resepctively. The 
theoretical en^ors correspond to the uncertainties due to the variation of (i) the Gegenbauer moments, the 
decay constants, (ii) the heavy-to-light form factors and the strange quark mass, and (iii) the wave function 
of the B meson characterized by the parameter Ag, the power corrections due to weak annihilation and hard 
spectator interactions described by the parameters Pa.h, ^a,h, respectively. To obtain the errors shown in 
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TABLE IV: CP-averaged branching fractions (in units of 10^^) and direct CP asymmetries (in %) of some 
selective B — > PP decays obtained in QCD factorization for three distinct cases: (i) without any power 
corrections, (ii) with power corrections to penguin annihilation, and (iii) with power corrections to both 
penguin annihilation and color-suppressed tree amplitudes. The parameters and 0a are taken from Table 
Hm pc = 13 and <^c = —70°. The theoretical en^ors con^espond to the uncertainties due to the variation of 
(i) Gegenbauer moments, decay constants, quark masses, form factors, the parameter for the B meson 
wave function, and (ii) Pa.h, <j>A.H, respectively. 



Mode y^/o pA,c,(t>A,c With Pa, (^A With Pa,c, 0a,c Expt. [3] 







,o 1+5.8+0.7 
^-'•^-3.5-0.7 


10 0+7.9+8.2 
^^•-'-4.8-6.2 


10 0+7.9+8.2 
^^••'-4.8-6.2 


19.4 ±0.6 






;r+2.8+0.3 

-'•-'-1.7-0.3 


o 4+3.8+3.8 
''•^-2.3-2.9 


Q ^+3.8+3.8 
''•"J-2.2-2.9 


9.8±0.6 


^{B- - 




M q+6.9+0.9 


91 7+9.2+9.0 
^^•'-6.0-6.9 


91 7+9.2+9.0 
^^•'-6.0-6.9 


23.1 ±1.0 


^{B- - 




Q 1+3.6+0.5 
^•^-2.3+0.5 


19 /r+4.7+4.8 
i^-O_3.0-3.7 


1 9 ;r+4.7+4.9 
^^•-'-3.0-3.8 


12.9±0.6 


^{B- - 




1 ^+1.1+0.3 
^•"-0.7-0.4 


9 4+1.8+1.3 
^•^-1.1-1.0 


9 4+1.8+1.3 
^•^-1.1-1.0 


2.36 ±0.27 






f. 9+0.4+0.2 
"•^-0.6-0.4 


7 A+0.4+0.7 
'•'-'-0.7-0.7 


7 A+0.4+0.7 
' •'^-0.7-0.7 


5.16±0.22 




> nPn^) 


A /1 9 +0.29+0. 18 
"•^^-0.11-0.08 


A c9+0.26+U.il 
'^■-'^-0.10-0.10 


1 1+1.0+0.7 
^•^-0.4-0.3 


1.55 ±0.19" 


m{B- - 




A Q+0.9+0.6 
^•^-0.5-0.3 


A Q+0.9+0.6 
^•^-0.5-0.3 


;r Q+2.2+1.4 
-'•^-1.1-1.1 


5 5Q+0.41 
-'•-'^-0.40 


^{B- - 




A 4+0.6+0.4 
^•^-0.3-0.2 


A r+0.6+0.5 
^•-'-0.3-0.3 


^ A+1.2+0.9 
-'•'^-0.6-0.7 


4.1±0.3 


Acpi^ ■ 


-^K-n+) 


4 A+0.6+1.1 
^•"-0.7-1.1 


7 4+1.7+4.3 
'•^-1.5-4.8 


7 4+1.7+4.3 
'•^-1.5-4.8 


-9.811:? 


Acp{^ - 




^•"-1.8-3.0 


A 7^+1.88+2.56 
U. '-J_o.94-3.32 


^^^•"-3.8-4.3 


-1±10 


Acp{B- 




A 79+0.06+0.05 
'^-0.05-0.05 


A 90+0.03+0.09 
'^•^°-0.03-0.10 


A 90+0.03+0.09 
'^•^°-0.03-0.10 


0.9 ±2.5 


Acp{B- 




7 + 1.6+2.3 
'•-'-1.2-2.7 


;r+1.3+4.9 
■'-'-1.8-4.6 


4 Q+3.9+4.4 
^•^-2.1-5.4 


5.0±2.5 


Acp{B- 




99 1+7.7+14.0 
^^•^-16.7- 7.3 


19 7+7.7+13.4 
^^•'-5.0-15.0 


1 1 n+ 8.4+14.9 
^^•"-21.6-10.1 


-37 ±9 


Acp{^ - 




9+0.4+2.0 
"•^-0.5-1.8 


17 A+1.3+4.3 
1 / .U_j 2-8.7 


17 A+1.3+4.3 


38±6 


Acp{^ - 


7r«7rO) 


4+ 6.8+34.8 
-'-'•^-10.6-37.7 


Of. Q+8.4+48.5 
^"•^-6.0-37.5 


;r7 9 + 14.8+30.3 
-"•^-20.8-34.6 


4^5+25 


Acp{B- 


^71-71°) 


"•"O-0.01-0.02 


'^•'^'J-0.01-0.02 


A 1 1+0.01+0.06 
'^•^ ^-0.01-0.03 


6±5 


Acp{B- 


^7r-T]) 


_i 1 4+1.1+2.3 
^^•^-1.0-2.7 


1 1 4+0.9+4.5 
^^•^-0.9-9.1 


^ A+2.4+ 8.4 
-'•'^-3.4-10.3 


-13±7 



"If an S factor is included, the average will become 1 .55 ± 0.35 . 



these tables, we first scan randomly the points in the allowed ranges of the above nine parameters (specifi- 
cally, the ranges p^-0.1 < Pa < Pa±O.I, (/>^-20° < ^a < ^(>a±20°, < p// < 1 and < 0// < 271 are used 
in this work, where the values of p^ and (j)^ are displayed in Table Ull)) and then add errors in quadrature. 
More specifically, the second error in the table is referred to the uncertainties caused by the variation of Pa^h 
and <j)A,H^ where all other uncertainties are lumped into the first eiTor. Power corrections beyond the heavy 
quark limit generally give the major theoretical uncertainties. 

A. Branching fractions 
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The B Kti decays are dominated by penguin contributions because of |V,„yjj*^| ^ |K-.vK*l 
and the large top quark mass. For the ratios defined by 



Rr 



R„ 



(4.1) 



T{B ^ K^n-) ' " " 2r(S0 ^ K^n^) ' 
we have 7?^ = ^« ~ 1 if the other quark-diagram amplitudes are negligible compared with P' . The current 
experimental measurements give Rc = l.\2± 0.07 and Rn = 0.99 ± 0.07. In QCDF we have = 1 . 15 ± 0.03 
and R„ = 1.12 it 0.03, which are consistent with experiment. 

From Table ITVl we see that the predicted rates for penguin-dominated B — > PP decays to the zeroth order 
of 1 /ni}, expansion are usually (30 ~ 45)% below measurements (see the second column of Table II VI). Also 
the direct CP asymmetry Acp{K^ 71^) is wrong in sign. We use penguin annihilation dictated by = 1.10 
and <^A = —50° to fix both problems. 



B^Kt] 



(') 



Among the 2-body B decays, B — > Kx]' has the largest branching fraction, of order 70 x 10 ^, while 
S§{B TqK) is only (1 - 3) x 10"^. This can be qualitatively understood as follows. Since the T\ — r\' 
mixing angle in the quark-flavor basis r]q = [uu + dd) / \/2 and T]^ = ss 

77 =cos0T]^ — sin(^T]i, T]' = sin0T]^ + cos^T]^,, (4.2) 

is extracted from the data to be = 39.3° 1251], it is clear that the interference between the B —>■ Kr\q 
amplitude induced by the b — > sqq penguin and the B —>■ Krjs amplitude induced by — > sss is constructive 
for B — > Krj' and destructive for B ^ rjK. This explains the large rate of the former and the suppression of 
the latter. However, most of the model calculations still fall short of the data for ^{B ^ Krj'). 

Many possible solutions to the puzzle for the abnormally large Krj' rate have been proposed in the past: 
(i) a significant flavor-singlet contribution il5l .l50ll . (ii) a lai^ge B —>^ rj' form factor 15 IH . (iii) a contribution 
from the charm content of the rj', (iv) an enhanced hadronic matrix element {0\sY5s\ri') due to the axial 
U(I) anomaly 15211. (y ) a large chiral scale Mq associated with the rjq |153|,|54|], (vi) a long-distance charming 
penguin in SCET 115 5n . and (vii) a large contribution from the two-gluon fusion mechanism 15611 . 

Numerically, Beneke and Neubert already obtained ^{B^ —>■ K^rj') ~ ^(50 X 10-*^) in QCDF using 
the default values = P// = l^. Here we found similar results 57 x 10^^ (53 x 10^^) with (without) 
the contributions from the "charm content" of the rj'. In the presence of penguin annihilation, we obtain 
^{B^ K^rj') ~ 78 X 10"'' (71 x 10"^) with (without) the "charm content" contributions. Therefore, the 
observed large B —>■ Krj' rates are naturally explained in QCDF without invoking, for example, flavor-singlet 
contributions. Data on S ^ Kri modes are also well accounted for by QCDF. 



B 



Tin 



From Table ITVl we see that power connections to the color-suppressed tree amplitude have almost no effect 
on the decay rates of penguin-dominated decays, but will enhance the color-suppressed tree dominated 
decay B — > n^^n^ substantially owing to the enhancement of |a2| ~ ^(0.6) [see Eq. (14.41 ) below]. Since 
\Pew/C\ is of order 0.06 before any power corrections, it is very unlikely that an enhancement of /few 
through New Physics effects can render c = C + Pew large and complex. Notice that the central values of 
the branching fractions of B° :/r"7r° measured by BaBar Q and Belle S, (1.83 ±0.21 ±0.13) x 10"^ 
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TABLE V: CP-averaged branching fractions (in units of of 5 ^ PP decays obtained in various ap 



proaches. The pQCD results are taken from 1161 
calculations for B Krj^'^ js 31.154 



62 



63 



63 



64j, |65|]. Note that there exist several pQCD 



6611 and here we cite the pQCD results with partial NLO corrections 



ll63h . There are two solution sets with SCET predictions for decays involving rj and/or rj' 115 5h . 



Mode 



QCDF (this work) 



pQCD 



SCET 



13.61'°? 
20.4+18^] 

7+6.0 



Expt. [3] 



B- 
B 

B 

B 

5° 



K-n+ 
K+K 



f9.2+9.0 
-6.0-6.9 
1 9 ^+4.7+4.9 
^'^•■'-S.O-S.S 
7.9+8.2 
4.8-6.2 
+3.8+3.8 
-2.2-2.9 
+2.2+1.4 



21.7^ 



19.3 



5.9 
7.0 
1.1 
1.8 
0.10 



9 1+1.0 



-1.1 
H0.7 
-0.7 
hO.7 



1.1- 
0.44 
0.7- 
+ 1.04 
-0.4-0.3 
+0.9+0.7 
-0.5-0.5 
+0.03+0.03 
-0.02-0.03 
0.8 
0.6 



8.71^:^ 



+6.7 



6-513.8 

0.29l°i° 
l.f 



0.046 
1.75 



20.8 ±7.9 ±0.6 ±0.7 
11.3±4.1±1.0±0.3 
20.1±7.4±1.3±0.6 
9.4 ±3.6 ±0.2 ±0.3 
5.2±1.6±2.1±0.6 
5.4±1.3±1.4±0.4 
0.84 ±0.29 ±0.30 ±0.19 
1.1±0.4±1.4±0.03 

1.0±0.4±1.4±0.03 



23.1 ±1.0 
12.9±0.6 
19.4 ±0.6 
9.8±0.6 

5 CQ+0.41 

5.16±0.22 
1.55±0.19 

1.36l°i? 



0.15 



-0.11 
-0.10 



96+°-2i 

"•^°-0.19 



B- 



B 



K X] 



K ri' 



50 ^ ^0^ 



50 ^ ^0^ 



B r] 

B Ti-r]' 

S° rjri 

5° ri'ri' 



9 + 1.8+1.3 
^■-'-1.1-1.0 



70 /I+61.2+26.4 
/o.t_26.8-19.5 



1.6 



-1.5+1.1 
-0.9-0.8 



■O"' lA 9+56.5+24.7 

'^•^-24.9-18.4 



^ 0+ 1-2+0.9 
■-'■'^-0.6-0.7 



o + 1.3+0.9 
■^•°-0.6-0.6 



n ^^+0.03+0. 13 
"•-"J-0.02-0.10 



n 49+0.21+0.18 
"•^^-0.09-0.12 



09+0.13+0.07 
"•-'"^-0.05-0.06 



n 0^+0.24+0.12 
"•-'°-0.10-0.08 



o + 1.2+2.7+1.1 
-'■^-0.9-1.2-1.0 



^1 n+13.5+11.2+4.2 
8.2- 6.2-3.5 



9 1+0.8+2.3+1.0 
^•^-0.6-1.0-0.9 



CA a+11.8+11. 1+4.5 
J\J.J_ 3^2- 6.2-2.7 



A ,+1.3+0.4+0.6 
^•^-0.9-0.3-0.5 



2.4l[Jj±0.2±0.3 



Q OO+0.04+0.04 I AC 
^•^^-0.03-0.03 =•= ^•^-' 



0.19 ±0.02 ±0.031°;°^ 



67+°-^2 

"•" -0.25 



0.18±0.11 



0.22 



-0.14+0.08 
-0.06-0.06 



0.11 



+0.12 

-0.09 



2.7 ±4.8 ±0.4 ±0.3 

2.3 ±4.5 ±0.4 ±0.3 
69.5 ±27.0 ±4.4 ±7.7 
69.3 ±26.0 ±7.1 ±6.3 

2.4 ±4.4 ±0.2 ±0.3 
2.3±4.4±0.2±0.5 

63.2 ±24.7 ±4.2 ±8.1 
62.2 ±23.7 ±5.5 ±7.2 
4.9±1.7±1.0±0.5 
5.0±1.7±1.2±0.4 
2.4 ±1.2 ±0.2 ±0.4 
2.8±1.2±0.3±0.3 

0.88 ±0.54 ±0.06 ±0.42 

0.68 ±0.46 ±0.03 ±0.41 
2.3±0.8±0.3±2.7 
1.3±0.5±0.1±0.3 

0.69±0.38±0.13±0.58 
1.0±0.4±0.3±1.4 
1.0±0.5±0.1±1.5 
2.2 ±0.7 ±0.6 ±5.4 

0.57 ±0.23 ±0.03 ±0.69 
1.2±0.4±0.3±3.7 



2.36 ±0.27 



71.1±2.6 



1 19+0.30 

^•^^-0.28 



66.1±3.1 



4.07 ±0.32 



2.7 



hO.5 
-0.4 



< 1.5 
1.2±0.4 

< 1.0 

< 1.2 

< 1.7 
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and (1.1 ±0.3 ±0.1) x 10^^ respectively, are somewhat different in their central values. The charged mode 
B 71^ tt" also gets an enhancement as its amplitude is proportional to c?] +a2- The prediction of QCDF or 
pQCD (see Table IVl) for — > tt+tt^) is slightly too lai^ge compared to the data. This is a long standing 
issue. One possibility for the remedy is that there exists nn nn meson annihilation contributions in 
which two initial quark pairs in the zero isospin configuration are destroyed and then created. Indeed, in 
the topolog 
flH.this 

model considered by Hou and Yang |60] and elaborated more by one of us (CKC) — > tt+tt^ and 
K^n^ rates are reduced and enhanced roughly by a factor of 2, respectively, through FSIs. It should be 
remarked that in the pQCD approach, it has been shown recently that the color-suppressed tree amplitude 
will be enhanced by a soft factor arising from the uncanceled soft divergences in the kj factorization for 
nonfactorizable hadronic B decays lioi ]. As a consequence, the B° n^n^^ rate can be enhanced to the right 
magnitude. 



the topological quai^k diagram approach, this corresponds to the vertical W-loop diagram 115911 . As shown in 
this additional long-distance contribution may lower the tt+tt^ rate. In the final-state rescattering 



The decays B — > K^K^ and B^ K^K^ receive b^d penguin contributions and B° — > K^K proceeds 
only through weak annihilation. Hence, the first two modes have branching fractions of order 10^^, while 
the last one is suppressed to the order of 10^^. 



TABLE VI: CP-averaged branching fractions (in units of 10 ^) of S ^ 7rT](') decays. 





n rj' 








BaBar 


3.5±0.6±0.2 [67] 


0.9±0.4±0.1 [68J 


4.00 ±0.40 ±0.24 [67] 


< 1.5 [68J 


Belle 


1.81°:^ ±0.1 M} 


2.8±1.0±0.3 |W] 


4.2 ±0.4 ±0.2 [70] 


<2.5 [TjJ 


Average 


o 7+0.5 
^•'-0.4 


1.2±0.4 


4.1 ±0.3 


< 1.5 



The decay amplitudes of B — > ttt] are 

V2A{B- ^ n-ri) « A^^^^ [dpu{a2 + M + + a^] ±A^^;, [5p„(ai ±j82)±a4^] , 
-2A(S° ^ Tr^T]) « A^^^ [V,(«2 - Pi)+2a^ + +A^^^„ [V(-a2 + ocP] , (4.3) 

with a4 = a4 + jSs and similar expressions for 5 — > ttt]'. It is clear that the decays 5" T]('^7r° have very 
small rates because of near cancelation of the the color-suppressed tree amplitudes, while the charged modes 
T]' 'tt^ receive color-allowed tree contributions. From the experimental data shown in Table IVTl it is clear 
that the BaBar's measurement of Si{B^ — > Tl^rj') » ^{B^ —>■ n^^rj') is in accordance with the theoretical 
expectation, whereas the Belle's results indicate the other way around. Nevertheless, BaBar and Bell agree 
with each other on ^{B tttj). QCDF predictions for B ttt]' ^ agree well with the BaBar data. As for 
the pQCD approach, it appears that its prediction for K^rj') is too small. At any rate, it is important 

to have more accurate measurements of B ^ Kri^'\ 
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TABLE VII: Same as Table IVl except for direct CP asymmetries (in %) of 5 — > PP decays obtained in 
various approaches. 



Mode 



QCDF (this work) 



pQCD 



SCET 



Expt. [3J 



B- 

B 

5° 

5° 

B 

5° 

5° 

B 



K-71+ 
K+K 



n 90+0.03+0.09 
"^•^"-o.os-o.io 

4 Q+3.9+4.4 



-7.4" 



-1.7+4.3 



+7 



'^-1.5-4.8 
_i A ^+2.7+5.6 
^'-'■"-3.8-4.3 
_A 1 1+0.01+0.06 
"•^^-0.01-0.03 
17 A+1.3+4.3 
^ '-"-1.2-8.7 
^7 9 + 14.8+30.3 
■"•^-20.8-34.6 



-6.4 



-0.8+1.8 
-0.6-1.8 



-10.0 





+0.7+1.0 

-0.7-1.9 



10 



18^20 

"-'-34 
II 

29 




<5 

-11±9±11±2 
-6±5±6±2 
5±4±4±1 
<4 

20±17±19±5 
-58±39±39±13 
I.l±0.4±1.4±0.03 

1.0±0.4±1.4±0.03 



0.9±2.5 
5.0±2.5 

-9.8+1:2 

-liblO 

6±5 

38±6 
40+25 

^-'-24 



12 



-17 

-18 



B- 



K r] 



B r]' 
B^^K^ri 

B^ ^n^rj 

B K^Tf]' 

^ k\' 



1 1 9+ 8.5+15.2 
-ii.z,_22.0-10.3 



A ^9+0.66+1. 14 
"•■'^-0.53-0.90 



-21.4 



8.6+11.8 
-22.9-11.3 



o A+0.6+0.7 
-'•"-0.5-0.8 



-5.0 



-2.4+ 8.4 
-3.4-10.3 



1 /;+5.0+ 9.4 
^■"-8.2-11.1 



9+2.8+24.6 
"-^ -^-5.0- 15.6 



_7 O + 1.0+17.6 
'•-'-1.8-14.0 



1 1 7+6.8+3.9+2.9 
'^^■'-9.6-4.2-5.6 



^ 9+1.2+1.3+1.3 
""•^-1.1-1.0-1.0 



19 7+4.1+3.2+3.2 
'^^•'-4.1-1.5-6.7 



9 O+0.5+0.3+0.2 
^••'-0.4-0.6-0.1 



07+8+4+0 
"-"-6-4-1 



OO+6+4+0 
"■-'■-'-4-6-2 



A0+ 9+3+1 
"^^-12-2-3 



0^+10+2+2 
"-'"- 9-1-3 



33±30±7±3 -37±9 
-33±39±10±4 

-10±6±7±5 1.3+};^ 
0.7±0.5±0.2±0.9 
21±20±4±3 
-18±22±6±4 
11±6±12±2 5±5 
-27±7±8±5 

5±19±21±5 -13±7 
37±19±21±5 
21±12±10±14 6±15 
2±10±4±15 
3±10±12±5 
-7±16±4±90 
-24±10±19±24 



5" rjri 



^ ri'ri' 



-63.5 



-10.4+ 9.8 

- 6.4-12.4 



^+7.2+3.8 
-^^•^-6.8-4.8 



-44 q+3-1+8.5 
^^■^-3.1-9.2 



-33 



+2.6+4.1+3.5 
-2.8-3.8-0.0 



77 4+0.0+ 6.9+8.0 
' '^^-5.6-11.2-9.0 



90 7+10.0+18.5+6.0 
^-'•'- 6.9-16.9-8.5 



-9±24±21±4 
48±22±20±13 

70±13±20±4 

60±11±22±29 
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B. Direct CP asymmetries 



For pc ~ 1-3 and 0c ~ —70°, we find that all the CP puzzles in B 
in fourth column of Table IIVI The corresponding a'2's are 



PP decays are resolved as shown 



02(7171) w 0.60e 



-(55° 



a2iK7i:)^0.5\e 



-;58° 



(4.4) 



They are consistent with the phenomenological determination of /T^ ) ~ a2/a\ from a global fit to the 
available data jlSh . Due to the interference between the penguin and the lai^ge complex color-suppressed 
tree amplitudes, it is clear from Table HVl that theoretical predictions for direct CP asymmetries now agree 
with experiment in signs even for those modes with the significance of Acp less than 3a. We shall discuss 
each case one by one. 



Acp{K-7r+) 

Neglecting electroweak penguin contributions, the decay amplitude of S" —>■ K 7r+ reads 

^ K-71+) = AMoci + oc^ + PD- 
Following l^, the CP asymmetry of S*' — > K^n^ can be expressed as 



Acp{b!^^K ;r+)/?FM = -2sin7lmrFM, 



(4.5) 



(4.6) 



with 



K-71^ 



FM 



''FM 



Am 



K^7l- 



1 - 2 COS /Re rpM + I ''fm | 



Ar 



ai{7iK) 



-al{7iK)-^'^{7iKy 



(4.7) 



where the small contribution from has been neglected and the decay ampUtude of B — > K'^tz^ is given in 
Eq. (14.1 lb . Theoretically, we obtain rpM = 0.14 for 7 = 67.8° with a small imaginary part and /?fm = 0.91, 
to be compared with the experimental value /?fm = 0.84 it 0.04. In the absence of penguin annihilation, 
direct CP violation of B^ K^7l^ is positive as Ima| ^ 0.013. When the power correction to penguin 
annihilation is turned on, we have Im(a| + jSg ) ^ —0.039 and hence a negative Acp{K^7l^). This also 
explains why CP asymmetries of penguin-dominated decays in the QCDF framework will often reverse 
their signs in the presence of penguin annihilation. 



Acp(^-7r°) 

The decay amplitude is 

V2A{B- ■ 



(4.8) 



If the color-suppressed tree and electroweak penguin amplitudes are negligible, it is obvious that the am- 
plitude of K^Ti^ will be the same as that of K^K^ except for a trivial factor of 1/V2. The CP asymmetry 
difference AAkk = Acp{K^ 71^) —Acp{K^7l^) arising from the interference between P' and C' and between 
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P^^ff and T' is expected to be small, while it is 0.148 ±0.028 experimentally iQ]. To identify the effect due 
to the color-suppressed tree amplitude, we write 

AAk^ = OmSt'^Zt-oml - 2sin7lmrc + • • • , (4.9) 

where the first term on the r.h.s. is due to the interference of the electroweak penguin with color-allowed 
tree and QCD penguin amplitudes and 

/III 



rc 



(■v) 



fKF^-{0)-al{7iK)-Pl{7iKy ^ • ^ 



The imaginary part of rc is rather small because of the cancelation of the phases between a2 and + j8|. 
When soft coiTcctions to a2 are included, we have rc ~ 0.078 — 0.063/ . It follows from Eq. ( I4.9I ) that 
will become of order 0.13 . 



As first emphasized by Lunghi and Soni ||72[1 . in the QCDF analysis of the quantity AAkh, although 
the theoretical uncertainties due to power corrections from penguin annihilation are large for individual 
asymmetries Acp{K^ Ti'^^) and Acp{K^ n'^), they essentially cancel out in their difference, rendering the 
theoretical prediction more rehable. We finA AAkk = [123^^11];^)%, while it is only (1.9lJJ-^+}-^)% in the 
absence of power corrections to ai or to the topological amplitude C . 

Acp{K^K^) andAcp{K°n-) 
The decay amplitudes are 

VlAiS^^K^Ti^) = A^^(-<-j8f)+A^^(5p„a2 + ^<Ew) = -/'' + c', 
AiB-^K^'n-) = A^^« + i3f)=y, (4.11) 

where the amplitudes p' = P' — ^^ew "'"^A' ^'^^ ~ ^' +^ew ^^'^^ been introduced in Sec. 1. CP viola- 
tion of B K^n^ is expected to be very small as it is a pure penguin process. Indeed, QCDF predicts 
Acp{K^^7l^) ~ 0.003. If c' is negligible compared to p', Acp(^*'7r") will be very small. Just as the previous 
case, the CP asymmetry difference of the K^^n^ and ^^71^ modes reads 

AA'j.^^AcpiK^K') -Acp{K^n-) = (0.57+°:«n:i4)% + 2sin7lmrc + • • • , (4.12) 

where the first term on the r.h.s. is due to the interference between the electroweak and QCD penguin 
amplitudes. To a good approximation, we have AA^jj ~ —^kk- This together with the measured value 
of AAjf;j and the smallness of Acp{K^Ti^) indicates that Acp{K^nP) should be roughly of order —0.15. 
Using Imrc ~ —0.063 as discussed before, it follows from the above equation Acp{K^ 'nP) is of order 
-11%. More precisely, we predict Acp(^" 71°) = (-10.6+^]!^:^)% and AA^^ = (-11.0t|^t|j)%, while 
they are of order 0.0075 and 0.0057, respectively, in the absence of pc and 0c. Therefore, an observation of 
Acp{K^n!^) at the level of— (10~ 15)% will be a strong support for the presence of power corrections to c' . 
This is essentially a model independent statement. 

Experimentally, the cunent world average —0.01 ±0.10 is consistent with no CP violation because the 



BaBar and Belle measurements -0.13 ±0.13 ± 0.03 IITJ and 0.14 ±0.13 ±0.06 IM, respectively, are of 



opposite sign. Nevertheless, there exist several model-independent determinations of this asymmetry: one 
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is the SU(3) relation Ar(7r*^7r*^) = —Ar{K^n^) i75h and the other is the approximate sum rule for CP rate 
asymmetries It^I 

Ar{K-n+) + Ar{K^n-)^2[Ar{K-7i^) + Ar{K^n'^)], (4.13) 

based on isospin symmetry, where Ar{Kn) = r{B Kn) — r(S —>■ Kn). This sum rule allows us to extract 
Acp(^"7r*') in terms of the other three asymmetries of K^Ti^ ,K^K^,K^7i^ modes that have been mea- 
sured. From the current data of branching fractions and CP asymmetries, the above SU(3) relation and 
CP-asymmetry sum rule lead to Acp{K^n^) = -0.0731°;°^^ and Acp{K^n^) = -0. 15 ± 0.04, respectively. 



An analysis based on the topological quark diagrams also yields a similar result —0.08 ~ —0.12 i77h . All 
these indicate that the direct CP violation Acp(^°7r") should be negative and has a magnitude of order 0.10 . 

The world average of Acp(B" -^K-T]) = -0.37 ±0.09 due to the measurements -0.36 ±0.11 ±0.03 



from BaBar 16711 and — 0.39±0.16±0.03 from Belle 117 Oil differs from zero by 4.1a deviations. The decay 



amplitude of S — > T] is given by III] 

VlAiB-^K-T^) = Ak^^ [5p„ a2 + 2a^] + V2A^^^ [5pup2 + oc^ + oc^ + P^] (4. 14) 

± VlA^r^^ [5pca2 + <] +A^^K [SpuioCl + J82) ± < ± p^] , 

where the flavor states of the T7 meson, qq = {uu + dd) / \/2, ss and cc are labeled by the T]^, and T]^, 
respectively. Since the two penguin processes b — > sss and b —>■ sqq contribute destructively to B ^ Kr] (i.e. 
A^^^ = X^^^'^'^ has an opposite sign to A^^^ and A^^^), the penguin amplitude is comparable in magnitude 
to the tree amplitude induced from b usu, contrary to the decay B — > Kr]' which is dominated by large 
penguin amplitudes. Consequently, a sizable direct CP asymmetry is expected in B K rj but not in 



The decay constants f^, and are given before in Eqs. (13.5! ) and (13.6b . Although —2 MeV 
is much smaller than f^'\ its effect is CKM enhanced by VcbV*s/{VubV*s). In the presence of penguin 
annihilation, Acp{K^ri) is found to be of order 0.127 (see Table HVl). When pc and (pc are turned on, 
Acp{K^ri) will be reduced to 0.004 if there is no intrinsic charm content of the 77. When the effect of 
is taken into account, Acp{K^ri) finally reaches the level of —11% and has a sign in agreement with 
experiment. Hence, CP violation in B K t] is the place where the charm content of the T] plays a role. 
Two remarks are in order. First, the pQCD prediction for Arp(K^ri) is very sensitive to niqq, the mass 



of the rjq, which is generally taken to be of order niji. It was found in II53I1 that for irigg = 0.14, 0.18 and 
0.22 GeV, Acp{K^ri) becomes 0.0562, 0.0588 and —0.3064, respectively. There are two issues here: (i) 
Is it reasonable to have a large value of iriqq ? and (ii) The fact that Acp{K^ri) is so sensitive to m^q 



implies that the pQCD prediction is not stable. Within the framework of pQCD, the authors of Il79ll rely 
on the NLO connections to get a negative CP asymmetry and avoid the aforementioned issues. At the 
lowest order, pQCD predicts Acp{K^ri) w 9.3%. Then NLO corrections will flip the sign and give rise 
to Acp{K^ri) = (—11.1^^^^^)%. In view of the sign change of Acp by NLO effects here, this indicates 
that pQCD calculations should be canied out systematically to NLO in order to have a reliable estimate of 
CP asymmetries. Second, while both QCDF and pQCD can manage to lead to a con^ect sign for Acp(^^T7), 
the predicted magnitude still falls short of the measurement —0.37 ±0.09. At first sight, it appears that the 
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QCDF prediction Acp{K^ri) = — 0.221+Q }g2 (see TablellVI) obtained in the leading 1 /m/, expansion ali^eady 
agrees well with the data. However, the agreement is just an accident. Recall that in the absence of power 
connections, the calculated CP asymmetries for K~n^ and Ti^n^ modes are wrong in signs. That is why 
it is important to consider the major power corrections step by step. The QCDF results in the heavy quark 
limit should not be considered as the final QCDF predictions to be compared with experiment. 

Acp{7i-ri) 

As for the decay B n^rj, it is interesting to see that penguin annihilation will flip the sign of 
Acp{7i^ri) into a wrong one without affecting its magnitude (see Table HVl). Again, soft corrections to 
02 will bring the CP asymmetry back to the right track. Contrary to the previous case of B — > K rj , the 
charm content of the rj here does not play a role as it does not get a CKM enhancement. 



Acp{n+n ) 



It is well known that based on SU(3) flavor symmetry, direct CP asymmetries in Kk and tik systems 
are related as 111] '■ 



-Ar(7r+7r 



Ar(^7r° 



-Ar(7r"7r' 



The first relation leads to Acp{k+k^) = [=^(^"7r+)/M7r+7r")]Acp(^"7r+) 
agreement with the current world average of 0.38 ± 0.06 ||3|]. 
The decay amplitude is 

A(S° ^ 71+71-) = A^^ [5p„(«i + iSi) + < + i3f + • • • ] , 



(4.15) 

0.37 , which is in good 
(4.16) 



which is very similar to the amplitude of the K 71+ mode (see Eq. ( 14.51 )) except for the CKM matrix ele- 
ments. Since the penguin contribution is small compared to the tree one, its CP asymmetry is approximately 
given by 



Acp{7i+7t ) Si 2 sin yim r;i 



with 



id) 



(4.17) 



(4.18) 



Numerically, we obtain Imr;i:;t= 0. 107 (—0.033) with (without) the annihilation term jSg . Hence, one needs 
penguin annihilation in order to have a correct sign for Acp(7r^7r^). However, the dynamical calculation of 
both QCDF and pQCD yields Acp (71+ 71") ss0.17 ~0.20. It is hard to push the CP asymmetry to the level of 
0.38 . Note that the central values of current B factory measurements of CP asymmetry: —0.25 ±0.08 ±0.02 
by BaBar H and -0.55 ± 0.08 ± 0.05 by Belle Q, differ by a factor of 2. 



Acp(7r"7r») 

Just like the K^r] mode, penguin annihilation will flip the sign of Acp(7r'^7r'') into a wrong one (see Table 
HVl) . If the amplitude c = C + Pew is large and complex, its interference with the QCD penguin will bring the 
sign of CP asymmetry into the right one. As mentioned before, IPew/C] is of order 0.06 before any power 



22 



corrections. It is thus very unhkely that an enhancement of Pew through New Physics can render c large and 
complex. For the 02(7171) given by Eq. (14.41 ). we find that Acp{7l'^7i^) is of order 0.55, to be compared with 
the current average, 0.43lo 24 iH]. 

Acp(7r~7r") 

It is generally believed that direct CP violation of B — > 71^ tt" is very small. This is because the isospin 
of the tt" state is / = 2 and hence it does not receive QCD penguin contributions and receives only the loop 
contributions from electroweak penguins. Since this decay is tree dominated, SM predicts an almost null 
CP asymmetry, of order 10^^ ~ 10^^. What will happen if 02 has a large magnitude and strong phase ? We 
find that power corrections to the color-suppressed tree amplitude will enhance Acp(7r^7r") substantially to 
the level of 2%. Similar conclusions were also obtained by the analysis based on the diagrammatic approach 



lll5h . However, one must be very cautious about this. The point is that power corrections will affect not only 
02, but also other parameters a, with / ^ 2. Since the isospin of tt^tt" is / = 2, the soft corrections to 02 
and Gi must be conspired in such a way that 71^71*' is always an / = 2 state. As explained below, there 
are two possible sources of power corrections to 02'. spectator scattering and final-state interactions. For 



final-state rescattering, it is found in i43ll that effects of FSIs on Acp(7r^7r*') are small, consistent with the 
requirement followed from the CPT theorem. In the specific residual scattering model considered by one of 
us (CKC) U, n can only rescatter into itself, and as a consequence, direct CP violation will not receive 
any contribution from residual final-state interactions. Likewise, if large pu and are turned on to mimic 
Eq. ^AA[ . we find Acp{7t^7t^) is at most of order 10"^. (The resuh of Acp{7i^7t^) in QCDF Usted in Tables 
ITVland rvIIl is obtained in this manner.) This is because spectator scattering will contribute to not only 02 but 
also a\ and the electroweak penguin parameters a-i^\Q. Therefore, a measurement of direct CP violation in 
B — > K^nP provides a nice test of the Standard Model and New Physics. 



CP asymmetries in pQCD and SCET 

For most of the B PP decays, pQCD predictions of CP asymmetries are similar to the QCDF ones 
at least in signs except for KK, K^7f , K~r]' , 7cri, rjrj' ,ri'ri' modes. Experimental measurements of Acp 
in 71^7], 71^7]' modes are in better agreement with QCDF than pQCD. It is known that power corrections 
such as penguin annihilation in QCDF are often plagued by the end-point divergence that in turn breaks the 
factorization theorem. In the pQCD approach, the endpoint singularity is cured by including the parton's 
transverse momentum. Due to a different treatment of endpoint divergences in penguin annihilation dia- 
grams, some of the CP puzzles do not occur in the approach of pQCD. For example, pQCD predicts the 
right sign of CP asymmetries for — > 7r*';r" and — > tt^tj without invoking soft corrections to 02- 

For decays involving T] and T]', there are two sets of SCET solutions as there exist two different sets of 
SCET parameters that minimize x^- It is clear from Table IVTIl that the predicted signs of CP asymmetries 
for K^7i^\ 7l^7l^\ 71^7] disagree with the data and hence the AA/fjr puzzle is not resolved. Also the predicted 
CP violation for K^n^ and K*^k^ is of opposite sign to QCDF and pQCD. This is not a surprise because the 
long-distance charming penguins in SCET mimic the penguin annihilation effects in QCDF. All the B-CP 
puzzles occurred in QCDF will also manifest in SCET. (The reader can compare the SCET results of Acp 
in Tables Ivn] (for B PP) and |XlIIl|XlV] (for B VP) with the QCDF predictions in the third column of 
Tables [TVl and 1X1 ) This means that one needs other power corrections to resolve the CP puzzles induced by 
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charming penguins. In tlie cun^ent phenomenological analysis of SCET |82], the ratio of ^ /T^ ^ is small 
and real to the leading order. This constraint should be released. 



C. Mixing-induced CP asymmetry 

Possible New Physics beyond the Standard Model is being intensively searched via the measurements of 
time-dependent CP asymmetries in neutral B meson decays into final CP eigenstates defined by 

where Am is the mass difference of the two neutral B eigenstates, Sf monitors mixing-induced CP asymme- 
try and A/ measures direct CP violation (note that C/ = —Acp). The CP- violating parameters C/ and Sf can 
be expressed as 

l-IA/P 21mA/- 

where 

In the standard model A/ w r]fe^^'^ for b ^ s penguin-dominated or pure penguin modes with f]f = 1 (— 1) 
for final CP-even (odd) states. Therefore, it is anticipated in the Standard Model that —rj/Sf ~ sin2j8 and 

The predictions of 5/ of B ^ PP decays in various approaches and the experimental measurements from 
BaBar and Belle are summarized in Table [Villi It is clear that rj'Ks appears theoretically very clean in 
QCDF and SCET and is close to sin2j8 = 0.672 it 0.023 determined from b — > ccs transitions jsj]. Note 
also that the experimental eiTors on Sr^i^s ^i'^ the smallest and its branching fraction is the lai^gest, making it 
especially suitable for faster experimental progress in the near future. 

Time-dependent CP violation in B° — > n'^Ks has received a great deal of attention. A correlation between 



^TfiKs ^^'^ Acpi^'^Ks) has been investigated in 18811 . Recently, it has been ai^gued that soft corrections to the 



color-suppressed tree amplitude will reduce the mixing-induced asymmetry ^^jo^^ to the level of 0.63 lllOll . 
However, we find that it is the other way around in our case. The asymmetry Sj^Of^^ is enhanced from 0.76 
to 0.791o'q41q o4 in the presence of power correction effects on a2- Our result of S^)]^^ is consistent with 
10. Si where power corrections were studied. ^ Although this deviates somewhat from the world average 
value of 0.57 ± 0.17 Hi, it does agree with the Belle measurement of 0.67 ± 0.3 1 ± 0.08 [85]. 

In sharp contrast to QCDF and SCET where the theoretical predictions for Sr^i^s ^i'^ ^^ry clean, the 
theoretical enws in pQCD predictions for both Sjjiks ^^id '^ijJfs arising from uncertainties in the CKM angles 
a and / are very lai^ge f;63i1 . This issue should be resolved. 

For the mixing-induced asymmetry in 5 — > 7r+;r^, we obtain S'^j+^j- = — 0.69lQ [QlQ Qg, in accordance 
with the world average of -0.65 ±0.07 U. For comparison, the SCET prediction —0.86 ± 0. 10 [550 is too 



^ Since power corrections will affect not only a2, but also other parameters a,- with ; ^ 2, we have examined such 
effects by using pn ~ 4.9 and 0// « —77° (see discussions after Eq. ( 13.201 ) ) and obtained the same result as before. 
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TABLE VIII: Mixing-induced CP violation Sf m B ^ PP decays predicted in various approaches. The 
pQCD results are taken from i6ll.l63l.l64l.l65n. For final states involving T] and/or T]', there are two solutions 
with SCET predictions i55n . The parameter T]/ = 1 except for Ks{'!i^ , rj , rj') modes where rjf = — I. Exper- 
imental results from BaBar (first entry) and Belle (second entry) ai^e listed whenever available. The input 
values of sin2j3 used at the time of theoretical calculations aie displayed. 



Mode 


QCDF (this work) 


pQCD 


SCET 


Expt. [83. 84. 85, 86, 87] 


Average 


sin2j3 


0.670 


0.685 


0.725 




ri'Ks 


n f;7+o.oi+o.oi 
"•"J'-o.oi-o.oi 


0.63^«|? 


0.706 ±0.008 
0.715 ±0.010 


0.57 ±0.08 ±0.02 
0.64 ±0.10 ±0.04 


U.jy ± u.u/ 


r]Ks 


7Q+0.04+0.08 

— U.UD — U.UO 


0.62;o|o 


0.69±0.16 
0.79 ±0.15 








n 7Q+0.06+0.04 
yj' 1 y fid. n OA 


0.74+002 


0.80 ±0.03 


0.55 ±0.20 ±0.03 
0.67 ±0.31 ±0.08 


A ^7 1 A 1 7 

U.J / ih U. 1 / 




n ^;q+0.08+0.19 

— U. lU — \).\jy 


-0.42+100 

—U.JO 


-0.86±0.10 


-0.68 ±0.10 ±0.03 
-0.61±0.10±0.04 


_A ^5 _|_ A 07 




n ns+o.06+0.19 

u.uo_Q j2_o.23 


067+° o*'5 


-0.90 ±0.24 
-0.67 ±0.82 








A if:+0.05+0.11 
"•^"-0.07-0.14 


067+0-°*"^ 
"•""'-0.011 


-0.96 ±0.12 
-0.60±1.31 








A 77+0.07+0.12 
"• ' ' -0.05-0.06 


A ror+0.004 


-0.98 ±0.11 
-0.78 ±0.31 








A 7^+0.07+0.06 
"• '°-0.05-0.03 


-0 131+0056 

"•^-"^-0.050 


-0.82 ±0.77 
-0.71 ±0.37 








A or +0.03+0.07 
"•"-'-0.02-0.06 


A qo+0.08 

"•^-'-0.12 


-0.59±1.10 
-0.78 ±0.31 







TABLE IX: Same as Table I VIII] except for A5/ for penguin-dominated modes. The QCDF results obtained 
by Beneke 19011 are displayed for comparison. 



Mode 


QCDF (this work) 


QCDF (Beneke) 


pQCD 


SCET 


Expt. 


Average 


With pc,0c 


W/o pc 


^'Ks 


oo+o °i 

"■""-0.01 


A A, +0.01 

"•"^-0.01 


A A, +0.01 

"•"^-0.01 


-0 06+"-^'' 

"•""-0.91 


-0.02 ±0.01 
-0.01 ±0.01 


-0.10±0.08 
-0.03 ±0.11 


-0.08 ±0.07 


r]Ks 


A ,2+0-09 
"•^^-0.08 


"•^^-0.03 


"•^"-0.07 


"•"'-0.92 


-0.04 ±0.16 
0.07 ±0.15 








12+°°^ 

"•^^-0.06 


09+° '" 

"•"^-0.06 


07+° °^ 

"•"'-0.04 


A AA+0.02 

0.06_()()3 


0.08 ±0.03 


-0.12±0.20 
0.00 ±0.32 


-0.10±0.17 



large and the theoretical uncertainty of the pQCD result — 0.42+H9 16111 is too large. For TrO^^C) modes. 
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SCET predictions are opposite to QCDF and pQCD in signs. For T]T]', the pQCD result is very small 
compared to QCDF and SCET. 

The reader may wonder why the QCDF result Sjj'Ks ~ 0.67 presented in this work is smaller than the 



previous result ^ 0.74 obtained in i89l . 



9C 



9111 . This is because the theoretical calculation of Sf depends on 
the input of the angle j8 or sin2j8. For example, sin2j3 0.725 was used in the earlier estimate of Sf around 
2005, while a smaller value of 0.670 is used in the present work. Therefore, it is more sensible to consider 
the difference 



ASf= -rifSf-sinip 



(4.22) 



for penguin-dominated decays. In the SM, Sf for these decays should be nearly the same as the value 
measured from the b — > ccs decays such as 5° — > J/yK^', there is a small deviation at most ^(0.1) 19311 . 
In Table IVIlTl we have hsted the values of sin2jS used in the theoretical calculations. Writing the decay 
amplitude in the form 



f) 



it is known that to the first order in rj = (A„A^ )/(AcAp 1941.19511 

ASf = 2|rj| cos2j8 sin/cosSj, 



(4.23) 



(4.24) 



with 8f = ai-g(A"/Ay^). Hence, the magnitude of the CP asymmetry difference A5j is governed by the size 
of A"/Ay^. In QCDF the dominant contributions to A" /A'^ are given by 19011 



A" 




i-n - [c] 


[-{al + r^al)]-[a"^R^,Ks] 


A^ 








A" 




[P"] + [c] [- 




A^ 




[PI 


[-{al + r^al)] 


A" 




[-P"] + [C] 




A^ 


7l»Ks 







(4.25) 



where R's are real and positive ratios of form factors and decay constants and we have followed 19011 to 
denote the complex quantities by square brackets if they have real positive parts. For rj'Kg, [—P] is enhanced 
because of the constructive interference of various penguin amplitudes. This together with the destructive 
interference between penguin and color-suppressed tree amplitudes implies the smallness of ASn'Ks- 
explained before, the penguin amplitude of rjKs is small because of the destructive interference of two 
penguin amplitudes [see Eq. (14.141 )1. This together with the fact that the color-suppressed tree amplitude 
contributes constructively to A"/A^ explains why ASj^Ks positive and sizable. 



^ The experimental value of sin 2/3 determined from all B-factory charmonium data is 0.672 ± 0.023 |l3|l. However, 
as pointed out by Lunghi and Soni one can use some observables to deduce the value of sin 2)3: CP -violating 
parameter Ek, AMs/AM^/ and Vcb from experiment along with the lattice hadronic matrix elements, namely, the kaon 
B-parameter Bk and the SU(3) breaking ratio ^j . A prediction sin2)3 = 0.87 ± 0.09 is yielded in the SM. If the ratio 
|Kzp/Ki| is also included as an input, one gets a smaller value 0.75 ±0.04. The deduced value of sin2j3 thus differs 
from the directly measured value at the 2a level. If the SM description of CP violation through the CKM-paradigm 
with a single Cf -odd phase is correct, then the deduced value of sin 2/3 should agree with the directly measured 
value of sin 2)3 in B-factory experiments. 
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Mixing-induced CP asymmetries in various approaches are listed in Table |IX] where the soft effects 
due to pc and 0c also displayed. In the QCDF approach, soft connections to the color-suppressed tree 
amplitude will enhance AS'^jO^^, shghdy from iff {0.09) to ^(0.12). It is clear that the OCDF results in the 
absence of power coiTcctions ai^e consistent with that obtained by Beneke ll90n . by us 11891] and by Buchalla 
et al. i9lll . For example, we obtained S-qiKg ~ 0.737 in 2005 and S-qix^ ~ 0.674 this time. But the value of 
^S^lKg remains the same as the value of sin2j8 has been changed since 2005. 



V. B^VP DECAYS 

Power corrections to ^2 for B —>■ VP and B —>■ VV are not the same as that for B — > PP as described by 
Eq. (14.41 ). From Table 1X1 we see that an enhancement of a2 is needed to improve the rates of B ^ p'^Tr" 
and the direct CP asymmetry of B° — > K*^ri. However, it is constrained by the measured rates of p^TT^ 
and p^tt" modes. The central values of their branching fractions are already saturated even for vanishing 
Pc{VP). This means that pc{VP) is prefeiTcd to be smaller than pc{PP) = 1.3 . In Table 1X1 we show the 
dependence of the branching fractions and CP asymmetries in B ^ VP decays with respect to Pa,c and 0a,c- 
The corresponding values of a2 for pc = 0.8 and 0c = —80° are 

a2{np)^0A0e-'^^\ a2(p?r) 0.38e"'^2° ^ 

a2{pK)K036e-'^'^\ a2(7r^*) w 0.39e"'^'° . (5.1) 

It is clear from Table 1X1 that in the heavy quark limit, the predicted rates for B — > ^*7r are too small by 
a factor of 2 3, while 3S{B —>■ Kp) are too small by (15 ^ 100)% compared with experiment. The rate 
deficit for penguin-dominated decays can be accounted by the subleading power corrections from penguin 
annihilation. Soft connections to a2 will enhance ^^{B p'^Tr") to the order of 1.3 x 10^^, while the BaBai^ 



and Belle results, (1.4±0.6±0.3) x 10"^ 1961] and (3.0±0.5 ±0.7) x 10"^ 1971] respectively, differ in their 



central values by a factor of 2. Improved measurements are certainly needed for this decay mode. 



A. Branching fractions 



B — > pn,a)K 

From Table Kill it is evident that the calculated B — > pn, con rates in QCDF are in good agreement with 
experiment. The previous QCDF predictions 111] for B — > pn (except B" — > Tl'^p*') are too large because 
of the large form factor Aq^{0) = 0.37 ±0.06 adopted in uJ]. In this work we use the updated sum rule 
result Aq''(0) = 0.303 ± 0.029 Q- It appeal's that there is no updated pQCD calculation for B ^ p;r and 
B (071. 

B^{p,(o,(^)ri(') 
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TABLE X: Same as Table HVl except for some selective B VP decays with pc = 0.8 and <pc = —80°. 



Mode 




W/o Pa,c,^a,c 


With Pa, 


With Pa,c, ^a,c 


Expt. [3] 






^ r+5.4+0.4 

—2.6—0.4 


^+5.7+7.4 
-2.8—4.5 


^+5.7+7.4 

—2.8—4.5 


8.6;?:? 






A 7+3.3+0.3 
— 1.7—0.3 


;r ^+3.5+4.3 
• —1.8—2.8 


^ 4+3.3+4.3 
— 1.7—2.8 


4.7 ±0.7 


m{B- - 




^ ^+6.1+0.7 
•^'^ 9 8 n ^ 

Z.o — V.J 


7 0+6.3+7.3 
•— 2.9— 4.4 


7 0+6.3+7.3 
'— 2.9— 4.4 


s.otJi 


^{B- - 


-K-p'') 


1 Q+2.5+0.3 
^ • in no 


o 0+2.6+2.9 

-^•-^ 11 17 

— i . i — I . / 


-r+2.9+2.9 
-'•'.'—I 2—1 8 


o 01 +0.48 

J • 1 r, / 


^(f - 


>K*-n+) 


T 7+0.5+0.4 

—V.J — U.'4 


Q 9+1.0+3.7 
y--^ in 

— 1 .U — J.J 


Q 9+1.0+3.7 

— i .U — J.J 


8.6+?-^ 

— 1 .U 


- 


> K*H^) 


1 ,+0.2+0.2 

• ^ —0.2—0.2 


^+0.4+1.7 

-^•"^ OS IS 

— U.J— I .J 


o ^+0.4+1.6 
-'■-J_0 4—1 4 


2.4 ±0.7 


m{B- - 




—0 9-0 6 


10.4+1-3+4.3 

— L .J — j.y 


10.4+I-3+4.3 

— i . J — J.y 




^{B- - 




9+0.4+0.3 
-'•'^-0 4-0 3 


f- 0+0.7+2.3 

—0.7—2.2 


f- 7+0.7+2.4 

' _Q 1 — 2 1 


6.9 ±2.3 


S8(t - 






1 5 4+V-7+9.4 
• —4 0—7 1 




15.9±1.0 


m(f - 


>coK^) 


Q+4.0+0.9 
— 1 6-0 4 


o Q+4.0+3.3 

— 1 6—2 2 


4 1+4.2+3.3 

— 1 7—2 2 


5.0 ±0.6 


m{B^ 


>p'^k'^) 


n 7f;+0.96+0.66 

—0 37-0 31 


n ro+O.88+0.60 
—0.32—0.22 


, + 1.7+1.2 

—0 6-0 6 


2.0 ±0.5'=' 


m{B- - 


^p-n') 


A A ^+1.2 + 0:9 

—0.9—0.5 


11 + 1.3+1.0 
—0.9—0.6 


11.8+I:«;i;^ 


10.9+}-^ 


m{B- - 


^p^Ti-) 


9 + 1.8+1.2 
—0.9—0.6 


^+1.8+1.2 
—0.9—0.6 


Q 7+2.7+1.7 
°- '-1.3-1.4 






>p^7l+) 


ic O+1.0+0.5 


i^r q+1. 1+0.9 
— 1.5—1.1 


ir Q+1. 1+0.9 
— 1.5—1.1 


15.7±1.8 


SS{B^ 


>p+n-) 


4+0.4+0.3 
° -0.7-0.5 


Q 9+0.4+0.5 
-0.7-0.7 


Q 9+0.4+0.5 
-0.7-0.7 


7.3±1.2 


Acpi^ ■ 


-^K-p+) 


1 O+0.7+3.8 
-0.3-3.8 


01 Q+11.5+19.6 
— 11.0—12.7 


01 Q+11.5+19.6 
— 11.0—12.7 


15±6 


Acp{^ - 




^ + 1.1+4.9 
— 1.2—4.9 


^ n+3-2+6.0 
—6.4—4.5 


7+1.2+8.7 
— 1.2—6.8 


_ 


Acp{B- 




r, 9/1+0.12+0.08 
n 1 S— 07 


n 97+0.19+0.46 
u.z, / _o.27— 0. 17 


n 97+0.19+0.46 
u.z, / —0.27-0.17 


-12±17 


Acp{B- 


^K-p') 


0+3.5 + 7.0 

—0 9—7 


^+16.1+30.0 
— 18 2—22 8 


4c 4+17.8+31.4 
— 19 4—23 2 


37±11 


Acp(f ■ 




1 c (,+0.9+4.5 

—0 7—4 7 


19 1+0.5+12.6 

U.J iO.U 


,9 ,+0.5+12.6 

—0 5— 16 


-23 ±8 


Acp(f ■ 




-12 0+2-^+ '^-^ 


n 07+1.71+6.04 

—0.89—6.79 


_, 7+1-8+9.1 
lu. / _2.8— 6.3 


-15±12 


Acp{B- 




n Q7+0.11+0.12 
—0 07-0 1 1 


aQ+0.04+0.10 
^•-''^— n n^- 1 7 


n aQ+0.04+0.10 
^'-''^—0 n^— 1 7 

u.VJJ VJ. IZ, 


-3.8±4.2 


Acp{B- 




,7 ;r+2.0+6.3 
— 1 3—8 


^ 7+0.7+11.8 
— 1.1 — 14.0 


1 zr+3. 1 + 11.1 
'— 1.7— 14.3 


4±29 


Acp(f ■ 




T 0+0-4+ 1-9 
—0 4—1 8 


n 90+0.51+2.00 
• —1 00—1 21 


;r+0.4+2.7 

OS 7 A 
W.J Z.M- 


19±5 


Acp(f - 




^ Q+1.9+3.4 
—2 3—4 1 


(. (:+4.7+6.0 

— J.'r — J.J 


A 7+1.8+5.5 

^* ' 1 S 8 

— i .0— J.O 


32 ±17 


Acp{B^ - 




9 0+2.4+9.9 
^•-'-3.7-9.2 


01 ^+13.3+21.5 
-'^■-'-12.5-30.9 


1 1 n+5-0+23.5 
^^•"-5.7-28.8 


-30 ±38 


Acp{B- 




^ 4+0.4+2.0 
-^•^-0.3-2.1 


1/; 3+1.1+ 7.1 
^"•-'-1.2-10.5 


Q 7+2.1+ 8.0 
^- '-3. 1-10.3 


2±11 


Acp{B- 


^p'k-) 


^ 7+0.5+3.5 
"•'-0.8-3.1 




Q 0+3.4+11.4 
^■''-2.6-10.2 


18^1^7 


Acp{B' - 


-^p-K+) 


o ^+0.2+1.0 
-'•■'-0.2-0.9 


4 4+0.3+5.8 
^•^-0.3-6.8 


A 4+0.3+5.8 
^•^-0.3-6.8 


11±6 


Acp(S° ■ 


-^p + K-) 


"•"-0.1-2.2 


99 7+0.9+8.2 
^^•'-1.1-4.4 


99 7+0.9+8.2 
^^•'-1.1-4.4 


-18±12 



"If an S factor is included, the average will become 2.0 ± 0.8. 
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The relevant decay amplitudes are 

V2A{B- ^p-ri) « Ap;,,, [5p„(a2 + i82)+2af + af] +A^„p [5p„(ai +p2) + «4] , 

V2A{B^ ^ 0T7) « A^^,0af + ^/iB^^^bl + VlB^^b'l (5.2) 

and similar expressions for T]'. It is clear that the decays B P^V^'^ have rates much larger than B^ — > 
p^rj^'^ as the former receive color-allowed tree contributions while the color-suppressed tree ampUtudes in 
the latter cancel each other. Both QCDF and pQCD lead to the pattern r{B ^ p^rj) > r{B p^T]'). 
This should be tested by more accurate measurements. The SCET prediction of ^{B^ ^ pi]') ~ 0.4 x 10^^ 
is far too small and clearly ruled out by experiment. Since the color-suppressed tree amplitudes in the 
decay B (orj'-'^ are added together, one should have r(B° G)T](')) > r(B° p*'t](')). It appears that 



SCET predictions for (p^jp", a))T]' il02|] are at odds with experiment. For example, solution I yields 
r(S° — > ft)T]') < r(S*' — > p^T]') in contradiction to the theoretical expectation and solution II gives r(S" — > 
p^rj ) > r{B P V ) disagreement with the data. 

The decays B° 01] ( ' are very suppressed as their amplitudes ai^e governed by ViibV*^i{a" — a'^)- For 
example, we obtain (prj) ^ 10^^ in the QCDF approach. Since the branching fraction of the corj 

mode is of order 10"^, it appears that the (j) meson can be produed from the decay B° (orj followed by 
(0 — (p mixing. This will be possible if (p is not a pure ss state and contains a tiny qq component. Neglecting 
isospin violation and the admixture with the p" meson, one can parametrize the ft)-0 mixing in terms of 
an angle 5 such that the physical ft) and ^ are related to the ideally mixed states co' = {uu + dd) / \f2 and 

= ss by 




cos 8 sin 8 
— sin 5 cos 5 




(5.3) 



and the mixing angle is about |5| ^ 3.3° il03h (see il04l] for the latest determination of 8). Therefore, the 



production of (^TJ through CO — <p mixing is expected to be 

^(5" ^ <pri)o,-^ iTiixing = ^{B^ (ori) sin^ 8^ 0.85 x 10"^ x (0.08)^ ~ 5.4 x 10"'^. (5.4) 

It turns out that the 0) — <p mixing effect dominates over the short-distance contribution. By the same token, 
the ft) — mixing effect should also manifest in the decay B —>■ <p7i^: 

S§{B- 07r")ffl-0 mixing = ^{B' (on-)^m^8K, 6.7 x 10"^ x (0.08)^ ~ 4.3 x 10"^ (5.5) 

For this decay, the short-distance contribution is only of order 2 x 10^^. 

B^K*K,KK* 

The decays B K*-K^,K*^K- and 5° K*^K^,K*^K^ are governed hy b ^ d penguin contribu- 
tions and 5° K*^K ,K* proceed only through weak annihilation. Hence, the last two modes are 
suppress ed rel ative to the first four decays by one order of magnitude. The recent preliminary measurement 



by Belle lll05h . ^{B- K*^K-) = (0.68±0.16±0.10) x 10"^ is in agreement with the QCDF prediction 



(see Table imi. 
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TABLE XI: Branching fractions (in units of 10 ^) of S — > VP decays induced by the b d (AS = 0) 



transition. We also cite the experimental data 1^13] and theoretical results given in pQCD |98, 
and in SCET mm 
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lod lioih 



Mode 


QCDF (this work) 


pQCD 


SCET 1 


SCET 2 


Expt. 


B- 




11.8+I^}1 

— 1 . i — i .4 


6^9 


Q+0.3+1.0 

— U. i — 1 .U 


11 4+0.6+1.1 

—U.O— U.y 


10.9^I-^ 


B 




Q 7+2.7+1.7 
— 1 . J — i .4 


5 ~6 


in 7+0.7+1.0 

i U. / (-17 AO 

— u. /—U.y 


J Q+0.2+0.8 
—0.1—0.8 


8.3+1-^ 

— 1.3 






9^ ,+1.5+1.4 

— 1.1— i .5 


18~45 


13.4+O-f }-2 

— U.D— 1.2 


16.8+°-^+j-^ 

— U.J— l.j 


23.0±2.3 


50 


P^Tl^ 


Q 9+0.4+0.5 
— u. /— u. / 




g Q+0.5+0.5 
' —0.5—0.5 


. ^+0.2+0.7 
— U. 1— u. / 


7.3±1.2 






15.9+11+?-? 

— 1 . J— 1 . i 




7 ^+0.3+0.8 
— U. 1 —U.o 


in 9+0.4+0.9 

lU.Z, /AC A Q 

—U.J— U.y 


15.7±1.8 


5« 


^ pOjrO 


1 + 1.7+1.2 
— U.o— U.O 


0.07 ~ 0.11 


9 r+0.2+0.2 

— U.l— U.i 


^ c+0.1+0.1 

—0.1—0.1 


2.0 ±0.5 


B 


con 


7+2.1+1.3 
— 1 0—1 1 


4~ 8 


7+0.4+0.7 

— U.J — U.U 


c+0.3+0.8 

—0 3-0 8 


6.9 ±0.5 




wtt" 


n 01 +0.02+0.04 
\j .yj L f\ r\r\ n n 1 
— u.uu — U.Ui 


0.10 ~ 0.28 


A AAAo+0.0299+0.0000 
\J.\j\J\Jj (\ (\(\[\r\ n nnAH 
— u . uuuu — u . uuuu 


A A1C+0.024+0.002 

U.UU f\ fiAn n nm 
— u.uuu — U.UUz 


<0.5 


B 


^ K*°K- 


n OA+O.20+0.31 
— U. 1 /— U.Zo 


0.32+012 

—U.U/ 


A /,Q+0.26+0.09 

— U.2U— U.Uo 


A ^ri +0.18+0.07 
U.J J- n 1 ^ A n*^ 
— U. Id— U.Uo 


0.68 ±0.19" 


B 


K*-K° 


4f,+0.37+0.42 

— U. i / — U.ZD 


0.211°}^ 


A r/, +0.26+0. 10 

— U.Zi — U.UO 


r^ ^1+0.21+0.08 

U.J 1 A 1 7 A 07 

— U. 1 / — U.U / 




5° 


K*+K- 


A Ao+0.01+0.02 

U.UO n n 1 n m 
— u.u 1 — u.uz 


— U . UD / 








5" 


K* K+ 


A A7+O.OI+O.O4 
— U.Ui — U.Uj 


0.017+««27 

— U . U i i 










K*^K° 


A 7A+0.18+0.28 
u. / u n 1 ^ n 

— U. i J— U.Zj 


0.24+00^, 

— u.uo 


A 4C+0.24+0.09 
yj.^u n 1 Q n n7 

— U. ly — U.U / 


A /I7+O.I7+O.O6 
U.T" / A I /I A A^ 

— U. 14— U.Uj 




5" 


K*^K° 


A /17+0.36+0.43 

U.T" / n 1 7 n 07 

— u. i / — u.z / 


0.49+0-^5 

— u.uv 


A ^1+0.24+0.09 

U.J 1 fi ofi fi fiQ 
— U.ZU — U.Uo 


A 40+0.20+0.07 
U.*+0 A 1 A A*i 
— U, 10 — U.UO 


< 1.9 


B 


—>■ <p7l^ 


w 0.043 ^ 


0.032+«°i2 


w 0.003 


« 0.003 


<0.24 


S« 




A A, +0.03+0.02 
U.UI a (w nm 

— U.Ui — U.Ui 


0.0068+[i|][J'« 

— U.UUUo 


« 0.001 


w 0.001 


<0.28 


B 


^p-T] 


o + 1.0+0.9 

— 0.0— u.y 




o Q+2.0+0.4 
— 1 . / — U.4 


3 3+1.9+0.3 

— 1.6—0.3 


6.9±1.0 


B 


^p-T]' 


c ^+0.9+0.8 
-^•"-0.5-0.7 


4.6^:1:^ 


A O7+2.46+0.08 
"•-"-0.22-0.07 


A /,/i+3. 18+0.06 
"•^^-0.20-0.05 


9 1+3.7 




^p^T] 


A ,A+0.02+0.04 
"■^"-0.01-0.03 




A A4+O.2O+O.OO 
"■"^-0.01-0.00 


A 14+0.33+0.01 
"•^^-0.13-0.01 


< 1.5 




^P^T]' 


A AQ+0. 10+0.07 
"•"^-0.04-0.03 


"•^"-0.05 


A 40+2.51+0.05 
"•^-'-0.12-0.05 


, A+3.5+0.1 
^•"-0.9-0.1 


< 1.3 


5« 


ft)T] 


A OC+O.65+0.40 
U. J _Q 26-0.24 


71+0.37 
"•'^-0.28 


A Q 1+0.66+0.09 
"•^^-0.49-0.09 


, 4+0.8+0.1 
^•^-0.6-0.1 


A Q4+0.36 
"•^^-0.31 




corj' 


A ;rQ+0.50+0.33 
"•■'^-0.20-0.18 




A ,0 + 1.31+0.04 
"•^°-0.10-0.03 


o 1+4.9+0.3 
-'•^-2.6-0.3 


^•"^-0.39 






w 0.005 ^ 


A All +0.062 
"•"^^-0.009 


K 0.0004 


w 0.0008 


<0.5 


50 




w 0.004 


017+"-'" 
"•"^ '-0.010 


0.0001 


« 0.0007 


<0.5 



"from the preliminary Belle measurement 1110511 . 
''due to the CO — mixing effect. 



B^K*7i,pK 

The relevant decay amplitudes are 

A{B^pK) = ApK{al-r^al + ^l + ---), 

A{B ^nK*) = A^K' (^4 ± 4* + ± • • • )• (5.6) 

Since the chiral factor is of order unity and r^' is small, it turns out numerically al{pK) ~ —al{nK*). 
Fortunately, ^^{pK) and ^^{nK*) are also of opposite sign so that penguin annihilation will contribute 
constructively. As noted before, in order to accommodate the data, penguin annihilation should enhance the 
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TABLE XII: Branching fractions (in units of 10 ^) of S — > VP decays induced by the b ^ s (AS = 1) t r ansi- 
tion. We also cite t he average of the experimental data and theoretical results given inpQCD |53, 106 1 

and inSCET mm . 



Mode 




QCDF (this work) 


pQCD 


SCET 1 


SCET 2 


Expt. 


B- 




^ 7+0.7+2.4 
-0.7-2.2 




A 9+2.2+0.8 
^•^-1.7-0.7 


^ ;r+1.9+0.7 
"•-'-1.7-0.7 


6.9±2.3 


B 




in 4+1.3+4.3 




^+4.7+1.7 
°-'-3.6-1.4 


q Q+3.5+1.3 
^■^-3.0-1.1 


Q Q+0.8 

^•^-0.9 


B^^ 




T -r+0.4+1.6 
■J -J -0.4- 1.4 




4 ^+2.3+0.9 
^•"-1.8-0.7 


T 7+1.4+0.5 
-"•'-1.2-0.5 


2.4 ±0.7 


B 


K*-n+ 


Q 9+1.0+3.7 
^•^-1.0-3.3 
o ^+2.9+2.9 
■-'••^-1.2-1.8 


6.oi:« 

5 1+^-1 


o 4+4.4+1.6 
°-^-3.4-1.3 
^ 7+2.7+1.0 
'J- ' -2.2-0.9 


Q ir+3.2+1.2 
^•-"-2.8-1.1 


8.61?:^ 

o o 1 +0.48 
-'•^^-0.46 


B 


p-K^ 


7 0+6.3+7.3 
'•°-2.9-4.4 


o 7+6.8 
»./_4.4 


Q 0+4.7+1.7 
^•-'-3.7-1.4 




8.0l}i 


B^^ 


p^K^ 


;r yl+3.4+4.3 
•^•^-1.7-2.8 




o jr+2.0+0.7 
-'■-'-1.5-0.6 


c Q+2. 1+0.8 
-"•O-l. 8-0.7 


4.7 ±0.7 




p+K 


o ^+5.7+7.4 
°-"-2.8-4.5 




Q 0+4.6+1.7 
^•°-3.7-1.4 


1 n 0+3.8+1.5 
^'-'•^-3.2-1.2 


8.61?:? 


B 

rS — > 


coK 

(OK 


A 0+4.4+3.5 
^•°-1.9-2.3 
A 1 +4.2+3.3 
^•^-1.7-2.2 


10.6ii°3^ 

Q Q+8.6 


1+2.4+0.9 
-'•^-1.9-0.8 
A 1+2.1+0.8 
^•^-1.7-0.7 


Q+2. 1+0.8 
-'•^-1.7-0.7 
A Q+1.9+0.7 
^•^-1.6-0.6 


6.7 ±0.5 
j.U ± u.u 


B 


(j)K- 


0+2.8+4.7 
0-0_2.7-3.6 
o 1+2.6+4.4 
°-^-2.5-3.3 


7.81^1 


Q 7+4.9+1.8 
^- -3.9-1.5 
Q ,+4.6+1.7 
^•^-3.6-1.4 


o ^+3.2+1.2 
O.U_2.7_i.o 

o ri+3.0+l.l 
^•U-2.5-1.0 


8.30 ±0.65 


B 


K* r] 


1 ;r 7+8.5+9.4 
iJ- '-4.3-7. 1 


00 1 O+0.26 


17 Q+5.5+3.5 
^ '-^-5.4-2.9 


IS A+4.5+2.5 


19.3 ±1.6 


B 


K*-ri' 


, 7+2.7+4.1 
^•'-0.4-1.6 


6.38 ±0.26 


A ;r+6.6+0.9 
^•-'-3.9-0.8 


A O+5.3+0.8 
^■^-3.7-0.6 


A Q+21 a 
^•^-1.9 


B^^ 




1 f,+7.9+9.4 
iJ-D-4.1-7.1 


-)-) T1+0.28 
ZZ. J l_o 29 


If, f:+5.1+3.2 
^0-0-5.0-2.7 


1^ 5+4.1+2.3 
^"•-'-4.3-2.0 


15.9±1.0 






, r+2.4+3.9 
^•■J-0.4-1.7 


o or+0.29 
-'•-'-'-0.27 


A 1+6.2+0.9 
^•^-3.6-0.7 


4 Q+4.7+0.7 
^•'J-3.4-0.6 


3.8 ±1.2^ 



"This is from the BaBar data \im] 
''This is from the BaBar data 110711 



Belle obtained an upper limit 2.9 x 10 ^ 
Belle obtained an upper limit 2.6 x 10^^ 



\m] 
\im 



rates by (15 ~ 100)% for pK modes and by a factor of 2 ~ 3 for K*n ones. A fit to the K*n and Kp data 
including CP asymmetries yields Pa{VP) « 1.07, (^a{VP) ^ -70°, Pa{PV) w 0.87 and MP^) ~ -30° as 
shown in Table Hill 

The pQCD predictions are too small for the branching fractions of K*^7l^ and K*^7l^, and too large for 
coK and coK^. 

B^^K 

A direct use of the parameter set pA{Py) ~ 0.87 and (j!>A(PV) -30° gives m{B- K-(j)) 13 x 10"^ 
which is too large compared to the measured value (8.30 ±0.65) x 10^^ This means that penguin- 
annihilation effects should be smaller for the (pK case. The values of Pa{K<P) and (j)A{K(j)) are shown in 
Table Hill It is interesting to notice that a smaller Pa for the ^ meson production also occurs again in VV 
decays. 

In the PP sector we leai^n that r(B Krj') » r{B Krj). It is the other way around in the VP sector, 
namely, r{B K*ri) » r{B K*r\'). This is due to an additional sign difference between a^{r\qK*) and 
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TABLE XIII: Same as Table IXj except for direct CP asymmetries involving b ^ d (AS' = 0) transitions. 



Mode 


QCDF (this work) 


pQCD 


SCET 1 


SCET 2 


Expt. 


B 




Q 7+2.1+ 8.0 
^- -3.1-10.3 


0~20 


1 ^ ;r+16.9+1.6 
^-'•-'-18.9-1.4 


19 -1+ 9.4+0.9 
^^•-'-10.0-1.1 


2±11 


B - 




Q 0+3.4+11.4 
^•°-2.6-10.2 


-20 ~0 


in + 13. 1+0.9 
lU.O_j27_o.7 


in 9 + 15.5+1.7 
^^•^-13.4-1.9 


18^1^7 






99 7+0.9+8.2 
'-1. 1-4.4 




n n+17.2+0.9 
^•^-16.7-0.7 


19 4+17.6+1.1 
-15.3-1.2 


-18±12 






A 4+0.3+5.8 
^-0.3-6.8 




11 0+17.5+1.2 
^^•°-20.0-l.l 


in 0+ 9.4+0.9 
^"•°-10.2-1.0 


11±6 






1 1 n+5-0+23.5 
^^•"-5.7-28.8 


-75 ~0 


_n g+21.4+0.1 
"•"-21.9-0.1 


o c+21.4+0.3 
-'•-'-20.3-0.3 


-30 ±38 


B -. 


con 


io 9+3.2+12.0 
lj.z_2.i_io.7 


~o 


n ^ + 19.1+0.1 

"•-'-19.6-0.0 


9 O + 13.4+0.2 
^•-'-13.2-0.2 


-4 ±6 




(071^ 


17 n+55.4+98.6 
1 / •W_22.8-82.3 


-20 ~ 75 


_n 4+24.0+1.1 
^•^-0.0-0.9 


on ^+ 79.1+3.4 
■''^-^-185.5-3.1 




B 


K*^K- 


Q+1. 1+2.8 
"•^-1.1-2.4 


^ Q+5.6+1.0+9.2+4.0 
"•^-5.3-0.3-6.5-6.0 


^+6.1+0.4 
-'•"-5.3-0.4 


_4 4+4.1+0.2 
^•^-4.1-0.2 




B 


K*-K^ 


7 0+5.9+ 4.1 
'•°-4.1-10.0 


f- ^+7.9+1.1+9.1+2.1 

"•■'-7. 3-1.4-7.7-3.9 


1 ^+2.6+0.1 
^•■'-2.3-0.1 


1 9 + 1.7+0.1 
^•^-1.7-0.1 






K*+K 


A 7+0.1+4.7 
^- ' -0.2-2.7 










5°^ 


K* 


c ^+0.2+7.0 
■'-'-0.2-5.5 














IT c+1.6+1.4 
^-'-'-1.7-2.3 




T ^+6.1+0.4 
■'•"-5.3-0.4 


_4 4+4.1+0.2 
^•^-4.1-0.2 








o ^+1.3+0.7 
-'•■'-1.7+2.0 




1 ;r+2.6+0.1 
^■^-2.3-0.1 


1 9 + 1.7+0.1 
^•^-1.7-0.1 




B 


(pTl^ 




















f- O+0.7+2.5 
"•■'-0.5-2.5 








B -. 


P T] 


r+0.4+6.5 
''•■'-0.4-5.3 


1 Q+0.1+0.2+0.1+0.6 
^•^-0.0-0.3-0.0-0.5 


_5 5+21.5+0.6 
"•"-21.3-0.7 


n 1+16.7+0.9 
^•^-15.8-0.8 


llibll 


B - 


p-T]' 


, +0.8+ 14.0 
^■^-2.2-11.7 


9c 1+0.8+2.1 
^-'•"-0.3-1.6-0.7-1.8 


10 0+66.5+2.8 
^^•"-37.5-3.1 


91 7+135.9+2.1 
^^•'- 24.3-1.7 


4±28 


B^^ 


pOrj 


o^: 9+3.7+10.4 
°"-^-5.8-21.4 


on ^+1.9+13.7+0.7+4.6 
"^•"-0.9- 3.9-0.1-9.0 


-46 7+170.4+2.9 
tu./_ 743_37 


00 0+66.9+3.1 
-'-'•-'-62.4-2.8 




B^^ 


pW 


^+4.5+39.5 
-'■'•■'-7.9-57.6 


nc 7+5.6+13.1+6.3+12.9 
'-'• '-4.8- 7.0-4.0- 9.9 


CI 7+103.3+3.4 
Ji-'_ 42.9-3.9 


r9 9 + 19.9+4.4 
-■^•^-80.6-4.1 




B^^ 


corj 


_44 7+13.1 + 17.7 
^•'- 9.9-11.6 


oo ^+1.0+0.8+5.9+3.9 
-'-'•■'-1.4-4.6-6.8-4.4 


n /I+30.7+0.9 
^•^-30.2-1.0 


n ^+17.8+0.9 
^•"-16.8-0.9 




B^^ 


corj' 


_41 4+2.5+19.5 
^^•^-2.4-14.4 


1^ n+0-1+3-3+2.2+1.7 
^"•"-0.9-3.9-3.2-2.0 


n+87.5+4.8 
^-'•"-38.8-5.1 


97 9 + 18.1+2.4 
z,/ .z,_29.7_2.2 





































0:4 (A'* T]^.) as discussed before. 

The QCDF prediction for the branching fraction of B ^ K*ri', of order 1.5 x 10^, is smaller compared 
to pQ CD and SCET. The experimental ave rages quoted in Table IXIII are dominated by the BaBar data 
lEoTn . Belle obtained only the upper bounds UM- ^{B^ ^ K*-ri') < 2.9 x 10"'' and ^{B- ^*°T]') < 
2.6 X 10^^. Therefore, although our predictions are smaller compared to BaBar, they are consistent with 
Belle. It will be of importance to measure them to discriminate between various model predictions. 



B. Direct CP asymmetries 



Acp{K*n) andAcp{Kp) 

First of all, CP violation for K*^n^ and p^K^ is expected to be very small as they are pure penguin 
processes (apart from a IV-annihilation contribution). From Table 1x1 we see that CP asymmetries for p^^K^, 
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TABLE XIV: Same as Table IXIII except for direct CP asymmetries (in %) involving A5 = 1 processes. 



Mode 



QCDF (this work) 



pQCD 



SCET 1 



SCET2 



Expt. 



B 




1 ^+3.1+11.1 

^•"-1.7-14.4 


_39+2l 


17 0+30.3+2.2 

^ '•°-24.6-2.0 


19 Q+12.0+0.8 
iz,.i/_j2.2-0.8 


4±29 


B 


K*^7i- 


n. 4+1.3+4.3 

^•^-1.6-3.9 










-3.8±4.2 


50 




in + 1-8+9.1 




c (^+^.5+0.5 
■'•"-8.4-0.5 


;r 4+4.8+0.4 
■'•^-5.1-0.5 


-15±12 


50 


-^K*-Ti+ 


19 1+0.5+12.6 
^^•^-0.5-16.0 


""-19 


11 9 + 19.0+1.3 
^^•^-16.2-1.3 


19 9 + 11.4+0.8 
^^•^-11.3-0.8 


-23 ±8 


B 




AC /I + 17.8+31.4 


7, +25 
'^-35 


Q 9 + 15.2+0.7 
^•^-16.1-0.7 


]f, 0+20.5+13 
1D.U_22.4-1.6 


37 ±11 


B 




r, O+0.2+0.5 
"•■-'-0.3-0.2 


1±1 








-12±17 


5° 




7+1.2+8.7 
"•'-1.2-6.8 


7+8 

'-5 


"•"-9.7-0.9 


o ^+4.8+0.3 
-'•-'-4.8-0.2 


6±20 


5° 


-^p+R- 


oi Q+11.5+19.6 
-'^•^-11.0-12.7 




7 1+11.2+0.7 
'•^-12.4-0.7 


Q ^+13.0+0.7 
^•"-13.5-0.9 


15±6 


B 


-^(OK 


99 ,+13.7+14.0 
^^•^-12.8-13.0 


39+15 

-'^-17 


, , ^+18.2+1.1 
ll.D_20.4-l.l 


19 + I6.6+O.8 
^^•-'-17.3-1.1 


2±5 


& 




A 7 + 1.8+5.5 
^•'-1.6-5.8 




9+8.0+0.6 
j.z_g2-0.6 


o 0+5.2+0.3 
-'•°-5.4-0.3 


32 ±17" 


B 




"•"-0.1-0.1 










-1±6 


5" 




Q Q+0.2+0.2 
"•^-0.1-0.1 


3+1 








23 ±15 


B- 




Q 7+3.9+6.2 
'-3.7-7. 1 


-24.57^0;^^ 


9 ^+5.4+0.3 
^•"-5.5-0.3 


, Q+3.4+0.1 
^•^-3.6-0.1 


2±6 


B 




^<r r + 10.1+34.2 
DJ.J_39 5_5Q2 


4.6011:1^ 


9 7+27.4+0.4 
^•'-19.5-0.3 


9 ^+26.7+0.2 
^•"-32.9-0.2 


-30tll 






o ;r+0.4+2.7 
-'•■^-0.5-2.4 


0.57 ±0.011 


, 1+2.3+0.1 
^•^-2.4-0.1 


n 7+1.2+0.1 
"• '-1. 3-0.0 


19±5 






^ + 10.7+33.2 
U.0_ 9 2-50.2 


-1.30±0.08 


Q (.+ 8.9+1.3 
^•"-11.0-1.2 


Q Q+6.2+0.9 
^•'^-4.3-0.9 


8±25 



''Note that the measurements of 52^11 ± 3 by BaBar II109I1 and 



9 ± 29 ± 6 by Belle lllOll are of opposite sign. 



p^K~ and K*^7l^ predicted in the heavy quark limit are all wrong in signs when confronted with exper- 
iment. For the last two modes, CP asymmetries are governed by the quantity rpM defined in Eq. (14.71 ) 
except that PP is replaced by VP or PV . Since a^{pK) and 6:4(71.^^*) are of opposite sign, this means that 
^Cp{P^K^) and Acp{K*^ 71^ ) should have different signs. This is indeed borne out by experiment (see Ta- 
ble lXIVl) . Numerically, we have a'^ipK) = 0.041 ±0.001/, a^ipK) = a^{pK) ± p^ipK) = 0.045 - 0.046/, 
a'^{7iK*) = -0.034± 0.009/ and a^{7iK*) = -0.066 ±0.013/. Therefore, one needs the ji^ terms (i.e. pen- 
guin annihilation) to get correct signs for CP violation of above-mentioned three modes. One can check 
from Eqs. (14.61 ) and (14.7b that Acp(p+.^^") is positive, while AcpiK*^^^) is negative. 
In order to see the effects of soft connections to a2, we consider the following quantities 

Acp{K*-n'^)-Acp{K*-n+)=Om6t''oZ-ofi 

Acp{K*''n'^)-Acp{K*''n-) = (-0.23iO:°}l°:U{)% ±2sin7lmrc(^*7r) ± • • • . 
defined in analog to AAkk and AA^jj with 

a2{K*7i) 



AAk'ji 



(5.7) 



rc{K*n) 



(■V) 



(■V) 



(5.8) 



fK'Ff''{0)-al{nK*)-^^,{7iK*] 

The first terms on the r.h.s. of Eq. ( 15.71 ) come from the interference between QCD and electroweak penguins. 
We will not consider similar quantities for Kp modes as the first term there will become large. In other 
words, as far as CP violation is concerned, /T*?! mimics Kn more than Kp. We obtain Im rc{K*n) = —0.057 
and Im rc(^p) = 0.023 and predict that AA^-*;: = (13.7;^^+^-^)% and AA^.^ = 



-ll.l+^-^+^-^)%, while it 
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is naively expected that K*^n^ and K*^n^ have similar CP-violating effects. It will be very important to 
measure CP asymmetries of these two modes to test our prediction. It is clear from Eqs. (15.71 ) and (14.61 ) (see 
also Tables ITVl and 1X1) that CP asymmetries of K*^n^ and K^^K^ are of order —0.10 and arise dominantly 
from soft connections to a2. As for Acp(^°p"), it is predicted to be 0.09 (ss —0.05) with (without) soft 
corrections to aj (cf. Table 1X1). 

Power coiTcctions to the color-suppressed tree amplitude is needed to improve the prediction for 
Acp(^*"t])- The current experimental measurement Acp(^*°T]) = 0.19 ±0.05 is in better agreement with 
QCDF than pQCD and SCET. 

In the pQCD approach, the predictions for some of the VP modes, e.g. AcpiK*-7i+), Acp(p°^") and 
Acp{p^K^) are very large, above 50%. This is because QCD penguin contributions in these modes are 
small, and direct CP violation arises from the interference between tree and annihilation diagrams. The 
strong phase comes mainly from the annihilation diagram in this approach. On the other hand, the predicted 
Acp{K*^ri) is too small. So fai- the pQCD results for Acp(^*T]^ ') are quoted from [53] where m,jy = 0.22 
GeV is used. Since the pQCD study of B — > Krj ''^ has been carried to the (partial) NLO and a drastic 
different prediction for Acp{K^ri) has been found, it will be crucial to generalize the NLO calculation to 
the K*ri^"> sector. 

We would like to point out the CP violation of (oK^. It is cleai^ from Table 1X1 that power coiTcction 
on a2 will flip the sign of Acp{(oK^) to a negative one. The pQCD estimate is similar- to the QCDF one . 
At first sight, it seems that QCDF and pQCD predictions are ruled o ut by the data Acp{(oK'^) = 0.032 it 



0.017. However, the BaBar and Belle measurements 0.52;°^2 _to.03 Ill09l] and -0.09 ±0.29 ±0.06 UM, 
respectively, are opposite in sign. Hence, we need to await more accurate experimental studies to test theory 
predictions. 

As for the approach of SCET, the predicted CP asymmetries for the neutral modes K*'^7i^,p^K^,coK^ 
and K*^ri have signs opposite to QCDF and pQCD. Especially, the predicted Acp{K*^ri) is already ruled 
out by experiment. 

Acp(p^) 

The decay amplitudes of B'^ ^ p'^K^ are given by 

A(B°^p-;r+) = A^p [5puOCi + + + ■ ■ ■] , 

A(B°^p+;r-) = Ap^[5puOCi + + + ■ ■ ■] . (5.9) 

Since the penguin contribution is small compared to the tree one, its CP asymmetry is approximately given 
by 



Acp{p 71+) 2sin7lmr;rp, 



with 



' Ttp 



A, 



{d) 



a'^{np)+l5^{np) 
cciinp) 



Acp{p^n ) PS 2sin7lmrp;r, 
a^{pK) + l5^{p7i) 



'pit 



Ac 



[d) 



OCiipTl) 



(5.10) 



(5.11) 



We obtain the values Imr^^p = 0.037 and Imrp;;; = —0.134. Therefore, CP asymmetries for p^n^ and p^7r+ 
are opposite in signs and the former is much bigger than the latter. We see from Table IXIlTl that the predicted 
signs for CP violation of p^Ti^ and p^Ti^ agree with experiment. The B —>■ p"?!^ decay amplitude reads 

A(B- ^ pOjr-) = Ap^ [5p„ai + < ± jSf ] +A^p [5p,a2 - < - p^] . (5. 12) 
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TABLE XV: Mixing-induced C P vi olation S/ in 5 — > VP decays predicted in vari ous a pproaches. The 
pQCD results are taken from i99lll06ll . There are two solutions with SCET predictions il02ll . The parameter 
T]y = 1 except for {<p,p,(o)Ks modes where rjf = —1. Experimental results from BaBar (first entry) and 
Belle (second entry) are listed whenever available. The input values of sin2j8 used at the time of theoretical 
calculations which are needed for the calculation of ASf are displayed. 



QCDF (this work) pQCD 



SCET 



Expt. IM>M>Ii3,I14>Ii5,I16,117,118J Average 



0.670 



0.687 



0.687 



n ^:qo +0.003+0.002 ATI -LA ni 

0-o92_o.ooo-o.oo2 0.71 ±0.01 



(^ 04^+0.05+0.04 

"•°^~0.05~0.06 



"■■^"-0.14-0.12 



0.84 



h0.03 
-0.07 



0.50 



-0.10 

-0.06 



0.51 
0.80 
0.85 
0.56 



0.69 
0.69 

0.05+0.02 
0.06-0.02 
0.02+0.01 
-0.02-0.01 
0.04+0.01 
-0.05-0.01 
0.02+0.01 
0.03-0.01 



0.26 ±0.26 ±0.03 



0.67 



-0.22 

-0.32 



0.55+g;|^±0.02 
0.11 ±0.46 ±0.07 
0.35l[J:^^± 0.06 ±0.03 
0.64;°;^^± 0.09 ±0.10 



Q 44+0.17 

"•^-0.18 



0.45 ±0.24 



54+0- 18 
'-'•■^^-0.21 



-0.24 



+0.15+0.20 
-0.14-0.22 



(^ 70+0.14+0.20 
'°-0.20-1.39 



0.51 



+0.08+0.19 
-0.07-0.32 



0.23 



h0.30 
-0.37 



A OA+0.04+0.24 
"•""-0.09-0.43 



-0 49+0-25 
"■^^-0.20 



_A if:+0.13+0.17 AC!Q+0.51 
"•^"-0.13-0.16 "-^^-0.66 



_n n +0.14+0.10 

"•^^-0.14-0.15 
_A 1Q+0.14+0.10 
"•^^-0.14-0.15 

_n 07+0.44+0.02 

"•°' -0.00-0.01 
r^ 79+0.36+0.07 
"•'^-1.54-0.11 
A o^:+0.15+0.03 
U.OD_2.03-0.07 
0.36+0.09 
0.44-0.15 
0.20+0.05 
1.73-0.09 
0.22+0.09 
-1.24-0.14 
0.19+0.10 
0.20-0.17 



0.29 
0.79 
0.38 
0.12 



A 1^+0.14+0.10 
""•^"-0.15-0.15 



(^ 90+O.I4+O.I6 
-U.zo_o.i3_o.i3 



0.77 



hO.22 
-0.53 



0.23] 



HO.59+0.10 
^-1.10-0.10 

n 97+0.17+0.09 
"•^'-0.33-0.14 



0.04 ±0.44 ±0.18 
0.17 ±0.57 ±0.35 



0.12±0.38 



As fai- as the sign is concerned, it suffices to keep terms in the first square bracket on the r.h.s. and obtain a 
negative Acp(p"7r^). By the same token, Acp{p^7i'^) is predicted to be positive. 

CP violation of B° p^n^ is predicted to be of order 0. 1 1 by QCDF and negative by pQCD and SCET. 
The current data are 0.10 ± 0.40 ± 0.53 by BaBar Imn and -0.49 ± 0.36 ± 0.28 by Belle wm . It is of 
interest to notice that QCDF and pQCD predictions for CP asymmetries of B ^ (p, ft))!]' ' are opposite in 
signs. 



C. Mixing-induced CP asymmetries 

Mixing-induced CP asymmetries Sf and ASf of B ^ VP in various approaches are listed in Tables IXVl 
and IXVII respectively. Just as the rj'Ks mode, <pKs is also theoretically very clean as it is a pure penguin 
process. Although the prediction of S,pKs ~ 0-69 has some deviation from the world average of 0.44+Qjg, 



it does agree with one of the B factory measurements, namely, 0.67^q32 by Belle (841] • In short, it appears 
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TABLE XVI: Same as Table IX Vl except for AS/ for penguin-dominated modes. The QCDF results obtained 
by Beneke 19011 are included for comparison. 



Decay 


QCDF (this work) 


QCDF (Beneke) 


pQCD 


SCET 


Expt. 


Average 


Withpc^c W/o pc 


^Ks 


022+°°°'^ 022+° °°^ 


02+" '" 


0.02 ±0.01 


~0 

~o 


-0.43 ±0.26 

o.o2i°:i 


-0 25+°+^ 


(oKs 
p'Ks 


n 17+0.06 n 1^+0.06 

"•^'-0.08 "■^-'-0.04 

-0 17+00'^ -0 11+007 


o.i3^«:«^ 
-0.081°:?^ 


15+003 
-0 19+0 10 


n 1 s+0.05 

11+002 

-0 13+*'°2 


-0.43 ±0.26 

-0 34+0-27 

"•-'^-0.31 

-0 05+''-2^ 


_n 14+0-26 

-0.10±0.17 



that the theoretical predictions of Sj for several penguin-dominated B PP,VP decays deviate from the 
world averages and hence may indicate some New Physics effects. However, if we look at the individual 
measurement made by BaBar or Belle, the theory prediction actually agrees with one of the measurements. 
Hence, in order to uncover New Physics effects through the time evolution of CP violation, we certainly 
need more accurate measurements of time-dependent CP violation and better theoretical estimates of S/. 
This poses a great challenge to both theorists and experimentalists. 



The ratio of A"/A'^ for the penguin-dominated decays {(}),(0,p'^)Ks has the expressions 19011 

[-K±^k)] 
[-K±'-k)]' 

[{al-r^al)] + [a^^R^Ks] 



A" 




i-P"] 




<t>Ks 




A" 




r']+[c] 


A^ 


(oKs 




A" 




[P"] - [C] 


A^ 




[PI 



(5.13) 



As discussed before, the quantity {a^ — r^a^) in above equations is positive and has a magnitude similar 
to |fl'4|. Since is larger than —a^^, ASj is positive for (oKs but negative for p^Ks and both have large 
magnitude due to the small denominator of ASf. From Table IXVll we see that ASo)Ks ~ ^(0-17), while 
^p°Ks ~ ^(—0.17). Effects of soft corrections on them are sizable. For example, ASpO^^ is shifted from 
f« —0.11 to « —0.17 in the presence of power corrections. This explains why our prediction of ASpOj^^ is 



substantially different from the Beneke's estimate ll90ll and our previous calculation [89]. 

For tree-dominated decays, so fai- there is only one measurement, namely, SpOj^o with a sign opposite to 
the theoretical predictions of QCDF and SCET. 

Time -dependent CP violation of the p+7r+ systems 

The study of CP violation for p^n^ and p^n^ becomes more complicated as p+7r+ are not 

CP eigenstates. The time-dependent CP asymmetries are given by 

T{B^{t) 



T{t{t) 



■p^n^] 



p^n^] 



r(s°(0 

{S ± AS) sin (Am? ) - (C ± AC) cos {Atnt ) , 



■p±7r+)±r(B0(?) ^p±;rT) 



(5.14) 
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TAB LE XVII: Various CP-violating param eters in th e decays S° p^TT^. SCET results ai^e quoted from 
1 102]. Experimental results are taken from 1 111 . 112 1 and the world average from Isj]. 



Parameter 


QCDF (this work) 


SCET 1 


SCET 2 


Expt. 


c 
s 

AC 
AS 


n 1 1 +0.00+0.07 
'-'•^^-0.00-0.05 
n nQ+0.00+0.05 
^•^^-0.00-0.07 

n r,/i+o.oi+o.io 

^•^^-0.01-0.09 
rv 9-^+0.02+0.02 
^•^"-0.02-0.02 
n 09+0.00+0.03 
^•^^^-0.00-0.02 


rv 19+0.04+0.04 
^•^^-0.05-0.03 

n ni +0.13+0.00 

^•^^^-0.12-0.00 
n 1 1 +0.07+0.08 
^•^^-0.08-0.13 

n 1 1 +0.12+0.01 

^■^^-0.13-0.01 
n /I7+0.08+0.05 
^•^'-0.06-0.04 


n 91+0.03+0.02 
^■^^-0.02-0.03 
rv rv 1+0.09+0.00 
^■^^-0.10-0.00 

n n 1+0 06+0.08 

^■^^-0.07-0.14 
rv 19+0.09+0.01 
"■^^-0.10-0.01 
r. /IO+0.05+0.03 
^■^•'-0.07-0.03 


-0.13 ±0.04 
0.01 ±0.07 
0.01 ±0.09 
0.37 ±0.08 

-0.04±0.10 



where Am is the mass difference of the two neutral eigenstates, S is referred to as mixing-induced 
CP asymmetry and C is the direct CP asymmetry, while AS and AC ai^e CP-conserving quantities. Defining 



= Ait 



and 



with qg/pg ~ e ^'^ , we have 



C + AC 



and 



p"7r+ 
— — — ) 

Pb ^+- 

|A+-p-|A+_p 

1 + |A+_|2 |A+_|2 + |A+_|2' 

2ImA+_ 



P 



A+_=A(F^^p+7r-), 



qg A. 



Pi, A_ 



1-|A+_| 



C-AC 



1 + |A_H 



|A-+P + |A_+| 



(5.15) 



(5.16) 



(5.17) 



S + ASi 
S-AS 



1 + 

2ImA_H 
1 + |A-+ 



2Im(f'2'/5A+_A;_; 

|A+_|2 + |A+_|2' 

2Im(e2'/3A_+Al+ 
|A_+|2 + |A_+|2 



(5.18) 



Hence we see that AS describes the sti'ong phase difference between the amplitudes contributing to B° 



p'^n^ and AC measures the asymmetry between r(S° p+Ti ) + r(S° ^ p 7r+) and r(5° ^ p 71+) + 



r(B°. 



P + TT 



Next consider the time- and flavor-integrated charge asymmetry 



|A+-|2 + |Ah 



- A- 



|A+_|2 + |A+_|2 + |A_+|2 + |A_ 



(5.19) 



Then, following 13411 one can transform the experimentally motivated CP parameters s^pji and Cpjt into the 
physically motivated choices 



Acp(p+7r ) = 



1 



+ 1 



with 



1b A- 



Pb Ah 



Acp{p 71+) 

Pb A-h 



+-|2 



1 



k+-|2 + i 



(5.20) 



(5.21) 
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Hence, 



r(f ^ p + n-) - r{B^ ^ p-n+) _ S^pn-Cpn-S^pnt^Cpn 
r(S°^p + 7r-)+r(SO^p-7r+) \-tS£pn-S^pnCpn 
r(f ^ p-7l+)-r{B^ ^ P + K-) _ £/pn + Cp^ + ^pnACpn 

r(f ^ p-n+)+r{B^ ^ p+n-) l+ACp„ + £/p„Cp„ 



(5.22) 



Therefore, direct CP asymmetries Acp{p^7l^) and Acp{p^7i^) are determined from the above two equa- 
tions and shown in TableslXland lXIIII Results for various CP-violating parameters in the decays B° p^TT^ 
are displayed in Table IXVlIl The CP- violating quantity ^/pj^ with the experimental value —0.13 ±0.04 is 
different from zero by 3.3a deviations. The QCDF prediction is in good agreement with experiment. 



VI. B^VV DECAYS 



A. Branching fractions 



In two-body decays S,,,/ PP,VP,VV, we have the pattern VV > PV > VP > PP for the branching 
fractions of tree-dominated modes and PP > PV ~ VV > VP for penguin-dominated ones, where B — > 
VP{PV) here means that the factorizable amplitude is given by {V{P)\J^\B){P(y)\J^\0). For example. 



p'n-) > 



(B 



71- Tt"), 



^{B- p-p") > m{B- p-K'^) > Sg{B- 

,^{B- ^ K^Ti-) > .^{B- ^ K*^7i-) ^ .^{B- ^ K*"" p-) > .^{B- ^K^'p-), (6.1) 

for tree- and penguin-dominated B decays, respectively. The first hierarchy is due to the pattern of decay 
constants /y > fp and the second hierarchy stems from the fact that the penguin amplitudes are proportional 
to + r^ag, ^4 + r^a6> fl^4 — + r^a^,, respectively, for B — > PP, PV, VP, VV. Recall that ~ 



1)» 

^ 



r^. There are a few exceptions to the above hierarchy patterns. For example, £^{B^ —>■ p^p*') < ^(5' 
n^n^) is observed. This is ascribed to the fact that the latter receives a lai^ge soft coiTcction to a2, while the 
former does not. 

There exist three QCDF calculations of B — > VV 112 ll . l22l . 12411 . However, only the longitudinal polarization 
states of B —>^ VV were considered in 1220. The analyses in 11211. 12411 differ mainly in (i) the values of 



the parameters p^ and 0a and (ii) the treatment of penguin annihilation contributions chai^acterized by the 
parameters jS, [see Eq. (12.101 )1 for penguin-dominated VV modes. Beneke, Rohrer and Yang (BRY) applied 
the values Pa{K*^) = 0.6 and (j)A{K*(j)) = -40° obtained from a fit to the data ofB ^ K*^ to study other 
B VV decays. However, as pointed out in 1240, the parameters Pa{K*p) 0.78 and ^a{K*p) ~ —43° 
fit to the data of B ^ K*p decays are slightly different from the ones Pa{K*(^) and (j)A{K*(j)). Indeed, we 
have noticed before that phenomenologically penguin annihilation should contribute less to <pK than pK 
and TtK*. This explains why the K*p branching fractions obtained by BRY ai^e systematically below the 
measurements. Second, as noticed in 12411 . there are sign errors in the expressions of the annihilation terms 



Aj'*^ and A'f obtained by BRY. As a consequence, BRY claimed (wrongly) that the longitudinal penguin 
annihilation amplitude jSj is strongly suppressed, while the {5^ term receives a sizable penguin annihilation 
contribution. This will affect the decay rates and longitudinal polarization fractions in some of B —>^ K*p 
modes, as discussed in details in 12411 . 
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In Table IX Vmi QCDF results are taken from 124(1 except that (i) a new channel coco is added, and 
(ii) branching fractions and fi for B ^ {p,K*)co decays are updated. ^ We see that the overall agreement 
between QCDF and experiment is excellent. In QCDF, the decay 5° — > wp" has a very small rate 

-2A(B" ^ cop'') « Ap„ [5p„{a2-pi)+2aP + a^] +A«p [5p„(-a2 -jSi) + <] , (6.2) 

due to a near cancelation of the the color-suppre ssed tree amplitudes. In view of this, it seems rather peculiar 
that the rate of p'^co predicted by pQCD il20(1 is larger than QCDF by a factor of 20 and exceeds the 
current experimental upper bound. Likewise, ^{B^ — > K*^p^) obtained by pQCD is slightly too large. 

We notice that the calculated p^p^ rate in QCDF is ^ = (0.88;^;^^+^;°^) x 10"'^ for 

pr = H, while BaBar and Belle obtained (0.92 ±0.32 ±0.14) x 10"^ and (0.4±0.4;U;2) x 10"^ 
II127I1 . respectively. Therefore, soft corrections to ai i.e. Pc(^V') should be very small for B^ — > p^p^. 
Consequently, a pattern follows: Power connections to a2 are large for PP modes, moderate for VP ones 
and very smaU for VV cases. This is consistent with the observation made in |0] that soft power correction 
dominance is much larger for PP than VP and VV final states. It has been argued that this has to do with 
the special nature of the pion which is a qq bound state on the one hand and a nearly massless Nambu- 
Goldstone boson on the other hand The two seemingly distinct pictures of the pion can be reconciled 
by considering a soft cloud of higher Fock states surrounding the bound valence quarks. From the FSI point 
of view, since B ^ P^P has a rate much larger than B ^ Ti^n^ , it is natural to expect that B —>■ n^n^' 
receives a large enhancement from the weak decay B ^ p^p^ followed by the rescattering of p^p^ to 
n^n'^ through the exchange of the p particle. Likewise, it is anticipated that B p^p" will receive a large 
enhancement via isospin final-state interactions from 5 —>^p^p^. The fact that the branching fraction of this 
mode is rather small and is consistent with the theory predic tion i mplies that the isospin phase difference of 
5q and ^2 and the final-state interaction must be negligible 112811 . 

Both ^*o^*o ^ 7i:*o^*- are b ^ d penguin-dominated decays, while B*' K* K*+ 

proceeds only through weak annihilation. Hence, their branching ratios are expected to be small, of order 
< 10"^. However, the predicted rates for and K*'^K*^ modes are slightly smaller than t he da ta. 

Note that a new Belle measurement of ^(S° k*^k*^) = (0.3 ±0.3 ±0.1) x 10"*^ < 0.8 x 10"^ JlOsIl is 
smaller than the BaBar result ^(B° f^*o^*o^ = (1.28lo 32) ^ ^^"^ U2M- Hence, the experimental issue 
with B — > K*K* decays needs to be resolved. 



B. Polarization fractions 

For charmless B — > VV decays, it is naively expected that the helicity amplitudes £/h (helicities h = 
0, — , ± ) for both tree- and penguin-dominated B VV decays respect the hierarchy pattern 

^,:^_:^^ = l:(^):(^)\ (6.3) 



Mb / \ nib 



Hence, they are dominated by the longitudinal polarization states and satisfy the scahng law, namely 112911 . 



/r-l-A = ^(4y ^ = l + ^f^), (6.4) 



m 



f\\ \mB 



^ The B ^ (0 transition form factors were mistakenly treated to be the same as that of B ^ p ones in the computer 
code of li^. Here we use the light-cone sum rule results from Isoll for Z? — > ft) form factors. 
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TABLE XVIII: CP-averaged branching fractions (in units of 10^^) and polarization fractions for B VV 
decays. For QCDF, the annihilation parameters are specified to be = 0.78 and <pA = —43° for K*p,K*K* 
and Pa = 0.65 and (pA = —53° for and K*(o by d efaii l t. The world a verages of experimental results 
are taken from iQ]. The pQCD results are taken from 1 119 . 



12C 



12lll23n . There are two distinct pQCD 
predictions for the branching fractions and longitudinal polarization fractions of B ^ K*{p,<p,(o) decays, 
depending on the type of wave functions. Numbers in parentheses are for asymptotic wave functions. Esti- 
mates of uncertainties are not available in many of pQCD predictions. 



Decay 










h 




QCDF 


pQCD 


Expt. 


QCDF 


pQCD 


Expt. 


B - 




^ „ i_A n 

^"■"-1.9-0.9 




^ , L 1 Q 

94 0+ 
-^^•"-2.0 


"•^"-0.01-0.02 




0.950 ±0.016 


t- 


-^p+p- 


l_l CiO A 

^-^•■^-2.6-1.5 


25.3+25-3 . 


24.2;^-i 


A QO+U.Ul+U.Ui 

^•^^-0.02-0.02 




0.978+«;°i 


t- 


-.pV° 


„ 11 C 1 1 1 

"•^-0.4-0.2 


„ 1 1 1 fl „ 

A q9+i.iu a 

"•^^-0.56 


"•'-"-0.28 


^•^^-0.04-0.37 


0.78 


"•'-"-0.15 


B - 


->■ p CO 


^ ^ n-i-l, 74-1 7 

^"•^-1.6-0.9 


19±2± 1 


15.9±2.1 


^•^"-0.01-0.03 


0.97 


0.90 ±0.06 


B - 


p^co 


r\ r\r^-\-0 07-1-0 '^ft 

"•"°-0.02-0.00 


1.9±0.2±0.2 


< 1.5 


v-^-l-0 1 1 -i-O 50 

"•-'^-0.25-0.36 


0.87 






-> coco 


-0.3-0.2 


1.2±0.2±0.2 


<4.0 


A Q^+U.U1+U.U4 
"■^ -0.01-0.20 


0.82 




B - 




A g+0.1+0.3 
"•"J-0.1-0.3 


A 40+0.12 

"•^°-0.08 


12 + 5 


A 4^+0.02+0.55 
"•^-'-0.04-0.38 


82 


75+016 
"• '■'-0.26 


t- 




A 1+0.0+0.1 
"•^-0.0-0.1 

A g+0.1+0.2 
"•"-0.1-0.3 


064+° °°^ 
35+°-" 


<2.0 

1 90+0.37 b 

^•^°-0.32 


1 

A jr9+0.04+0.48 
"•■'^-0.07-0.48 


0.99 
0.78 


o.so^S:!^ 


B - 


-^K*^p- ' 


Q 9 + 1.2+3.6 
'^•^-l. 1-5.4 


17 (13) 


9.2±1.5 


A 40+0.03+0.52 
"•^°-0.04-0.40 


0.82 (0.76) 


0.48 ±0.08 


B - 


-^K*-p^ 


c ^+0.6+1.3 
■^•■'-0.5-2.5 


9.0 (6.4) 


<6.1 


A ^7+0.02+0.31 
-0.03-0.48 


0.85 (0.78) 


A 95+0.06 d 

"•^"-0.16 




^K*-p+ 


o Q+1. 1+4.8 
°-^-1.0-5.5 


13 (9.8) 


< 12 


A ro+0.02+0.45 
"•-^-'-0.03-0.32 


0.78 (0.71) 




^- 




A ft+0.6+3.5 
^•"-0.5-3.5 


5.9 (4.7) 


3.4±1.0 


A oQ+0.00+0.60 
"•-'^-0.00-0.31 


0.74 (0.68) 


0.57 ±0.12 


B - 




iA A+1.4+12.3 
^"•"-1.3- 6.1 


/ 


lO.Oibl.l 


A 4Q+0.04+0.51 
"•^^-0.07-0.42 


/ 


0.50 ±0.05 


5°- 




Q ir+1.3+11.9 
^•■'-1.2- 5.9 


/ 


9.8±0.7 


A cA+0.04+0.51 
"•-'"-0.06-0.43 


/ 


0.480 ±0.030 


B - 




o A+0.4+2.5 
-0.3-1.5 


7.9 (5.5) 


<7.4 


A ^7+0.03+0.32 
"•"'-0.04-0.39 


0.81 (0.73) 


0.41 ±0.19 


5°- 


-^K*°co 


9 c+0.4+2.5 
^•J-0.4-1.5 


9.6 (6.6) 


2.0±0.5 


A ro+0.07+0.43 
"•■^ ''-0.10-0.14 


0.82 (0.74) 


0.70 ±0.13 



"There exist several pQCD calculations for pp modes il2(All2llll23[l . Here we cite the NLO resul 



pmlllll. 



''This is from the BaBar data lll24ll . The Belle's new measurement yields (0.3 ±0.3 ±0.1) x lO'^ fllOSH . 
'^^This mode is employed as an input for extracting the param eters and 0^ for B K*p decays. 
■^A recent BaBar measurement gives fi{K*^ p^) = 0.9 ± 0.2 I125j], but it has only 2.5c7 significance. 



■^This mode is employed as an input for extracting the parameters p^ and for B 
^See footnote 6 in Sec. VLB. 



■K*<p decays. 



with fL,f±,f\\ and being the longitudinal, perpendicular, parallel and transverse polarization fractions, 
respectively, defined as 



|2 



r |4)|2 + |^-|2 + |^^|2' ^6-^) 

with a = L, ||,±. In sharp contrast to the pp case, the large fraction of transverse polarization of order 
0.5 observed m B ^ K*<p and B —>^ K*p decays at B factories is thus a surprise and poses an interesting 
challenge for any theoretical interpretation. Therefore, in order to obtain a large transverse polarization in 
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B K*^,K*p, this scaling law must be circumvente d in o ne way or another. Various mec hanisms such 
as siza ble p enguin-induced a nnihilation con t ribut ions 112911 . final-state interactions i43l. 1 13011 . form-factor 
tuning ilBin and new physics 11132 , 
puzzle. 
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134j. 113511 have been proposed for solving the B — > VV polarization 



As pointed out by Yang and one of us (HYC) ll24ll . in the presence of NLO nonfactorizable corrections 
e.g. vertex, penguin and hai^d spectator scattering contributions, effective Wilson coefficients a'l are helicity 
dependent. Although the factorizable helicity amplitudes X^, X and X+ defined by Eq. (12.41 ) respect the 
scaling law (16.31 ) with Aqcd/'m/jj replaced by 2my /niB for the light vector meson production, one needs to 
consider the effects of helicity-dependent Wilson coefficients: .s^^ j s^^ = /(ar)X^/[/(a°)X°]. For some 
penguin-dominated modes, the constructive (destructive) interference in the negative-heUcity (longitudinal- 
hehcity) amphtude of the B VV decay will render f{a7) ^ so that is comparable to iz/" 

and the transverse polarization is enhanced. For example, fL{K*^p^) '-^ 0.91 is predicted in the absence of 
NLO corrections. When NLO effects ai^e turned on, their corrections on will render the negative helicity 
amphtude (S" K*^p'^) comparable to the longitudinal one K*^p^) so that even at the short- 

distance level, fi for B — > K*^p^ can be as low as 50%. However, this does not mean that the polarization 
anomaly is resolved. This is because the calculations based on naive factorization often predict too small 
rates for penguin-dominated B — > VV decays, e.g. B —>■ K*<p and B — > K*p, by a factor of 2 ~ 3. Obviously, 
it does not make sense to compare theory with experiment for fij as the definition of polarization fractions 
depends on the partial rate and hence the prediction can be easily off by a factor of 2 3. Thus, the first 
impo rtant task is to have some mechanism to bring up the rates. While the QCD factorization and pQCD 
1 13611 approaches rely on penguin annihilation, soft-collinear effective theory invokes charming penguin 



1711 and the final-state interaction model considers final-state rescattering of intermediate charm states 1431 
130l 113711 . A nice feature of the {S — P){S + P) penguin annihilation is that it contributes to and £/~ 
with similar amount. This together with the NLO corrections will lead to fi ~ 0.5 for penguin-dominated 
VV modes. Hence, within the framework of QCDF we shall assume weak annihilation to account for 
the discrepancy between theory and experiment, and fit the existing data of branching fractions and fi 
simultaneously by adjusting the parameters and <pA- 

For the longitudinal fractions inS K*p decays, we have the pattern (see also 1I21I]) 

fLiK*-p'') > MK*-p+) > fL{K*'p-) > fL{K*^p''). (6.6) 

Note that the quoted experimental value fL{K*-p^) = 0.96;o'|'g in Table IXVlII was obtained by BaBar in 
a previous measurement where ic'*^p" and 7c'*^/o(980) were not separated lll38ll . This has been overcome 
in a r ecent BaBar measurement, but the resultant value /i(^*^p") = 0.9 ±0.2 has only 2.5a significance 
1 12511 . At any rate, it would be important to have a refined measurement of the longitudinal polarization 



fraction for K*^p and K p and a new measurement of fL{K*^p^) to test the hierarchy pattern (16.61 ). 
In the QCDF approach, we expect that the b d penguin-dominated modes K*^^K*^ and K*^K*^ have 
~ 1/2 similar to the AS = 1 penguin-dominated channels. However, the data seem to prefer to fi ~ 

1^(0.75 — 0.80). Due to the near cancelation of the color-suppressed tree amplitudes, the decay B^ — > cop^ 

is actually dominated by b ^ d penguin transitions. Hence, it is expected that fL{p^co) ~ 0.52. It will be 

interesting to measure ft for this mode. 
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TABLE XIX: Direct CP asymmetries (in %) ofB^VV decays. The pQCD results are taken fro m j 12011 
for pp,p<j),K*K*. CP asymmetries of K*p and K*(0 in the pQCD approach are shown in Fig. 4 of ill9(1 as 
a function of y and only the signs of Acp{K*p) and Acp{K*(o) are displayed here. Note that the definition 
of Acp in ill9(1 has a sign opposite to the usual convention. 



Decay 



QCDF (this work) 



pQCD 



Expt. [3J 



B- 




0.06 





-5.1 ±5.4 




-^p+p- 


A+O+3 
^-0-3 


-7 


6±13 






-"-'-16-26 


80 






p CO 


8+1+3 


-23 ±7 


-20 ±9 




-^p'^co 


0+2+51 








(0(0 


on+15+16 

■"•"-M-IS 






B 




1^+1+17 


-15 




t 


K* K*+ 





-65 




—0 




. 1 1 1 A 

-14^}+^ 







B- 


^K*'^p- 


-0 3+°+2 


+ 


-1±16 


B 


^K*-p^ 


4^+6+12 
^•'-3-28 


+ 


2oii 




-^K*-p + 


00 + 1+ 2 
-'^-3-14 


+ 




t 




1 jr+4+16 
iJ-8-14 




9±19 


B^ 




0.05 




-1±8 






n O+0+0.4 




1±5 


B 


^K* (0 


56+3+ 4 


+ 


29 ±35 




-^K*°(o 


0+9+ 5 
^-'-5-18 


+ 


45 ±25 



For AS' = 1 penguin-dominated modes, the pQCD approach predicts ft ~ 0.70 — 0.80 . 



C. Direct CP asymmetries 



Direct CP asymmetries of S ^ VV decays are displayed in Table IXIXl They are small for color-allowed 
tree-dominated processes and large for penguin-dominated decays. Direct CP violation is very small for the 
pure penguin processes K*^p^ and K*(l). 



Early pQCD calcul ations of B ^ K*<j) tend to give a large branching fraction of order 15x10 and the polarization 



fracti on fi ^ 0.75 ll36l 113911 . Two possible remedies have been considered: a s mall f orm factor Aq*^ (0) = 0.32 
1 13111 and a proper choice of the hard scale A in B decays lll2lll . As shown in 112111 . the branching fraction of 
becomes 8.9 X lO"*^ and /l - 0.63 for A 1.3 GeV. 
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D. Time-dependent CP violation 



In principle, one can study time-dependent CP asymmetries for each helicity component, 

m)-r(B°(0^m) 



T{B\t) ^ nn) +r(so(o ^ nn; 

5/, sin ( Amf ) — C/, cos ( Amf ) . 



(6.7) 



Time-dependent CP violation h as been rn easured for the longitudinally polarized components of 
p^p^ and p^p" with the results Judluill : 



= -0.05 ±0.17, 
-0.3±0.7±0.2, 



c 



= -0.06±0.13, 
0.2 ±0.8 ±0.3. 



In the QCDF approach we obtain 

^{p+p-)L = {24Jtiitii)xlO-\ 

(0-6-0.3-0.3^ 



;+1.3+0.8^ ^ 



cP^P" 



n 1Q+0.01+0.09 rP 

"■^^-0.00-0.10' "-L 



. A 1 A+0 05+0.50 r 
^•^"-0.11-0.48' 



(6.8) 



n 1 1+0.01+0.11 

— U.11_QQJ_Q04, 

-0 53 +0-23+0.12 eg Qx 



As pointed out in IM, since (see Eq. (33) of UM and Eq. (106) of IE) 



sin 2a ± 2rp cos 5p sin 7 cos 2a ± ^(r 



(6.10) 



with P =\T\rp cos 5p and a = n — j5 — y, the measurement of can be used to fix the angle 7 with good 
accuracy. For the QCDF predictions in Eq. dm we have used j8 = (21.6lo^)° 7 = (67.8t3 9)° ibj- 



VII. CONCLUSION AND DISCUSSION 

We have re-examined the branching fractions and CP-violating asymmetries of charmless B — > 
PP, VP, VV decays in the framework of QCD factorization. We have included subleading I /mi, power 
corrections to the penguin annihilation topology and to color-suppressed tree amplitudes that are crucial 
for explaining the decay rates of penguin-dominated decays, color-suppressed tree-dominated n'^n^\ p'^Ti" 
modes and the measured CP asymmetries in the Bu.d sectors. A solution to the AAkk puzzle requires a large 
complex color-suppressed tree amplitude and/or a large complex electroweak penguin. These two possibil- 
ities can be discriminated in tree-dominated B decays. The CP puzzles with n^rj, K^n^ and the rate deficit 
problems with n^n^, p^n^ can only be resolved by having a large complex color-suppressed tree topology 
C. While the New Physics solution to the B Kit CP puzzle is interesting, it is irrelevant for tree-dominated 
decays. 

The main results of the present paper are: 
Branching fractions 

i). The observed abnormally large rates of S — > Kr]' decays are naturally explained in QCDF without 
invoking additional contributions, such as flavor-singlet terms. It is important to have more accurate 
measurements of B ^ ;rT](') to confirm the pattern 3S{B^ ^^V') ^ ^(5" 7r"T]'). 
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ii) . The observed large rates of the color-suppressed tree-dominated decays — > 7l^7l'^,p^n^ can be 

accommodated due to the enhancement of \a2{7i7i)\ ~ ff{0.6) and \a2{7cp)\ ~ ff{OA). 

iii) . The decays — > 0T] and B (pn^ are dominated by the co — <p mixing effect. They proceed 

through the weak decays (orj and B — > (on , respectively, followed by ft) — mixing. 

iv) . QCDF predictions for charmless B —>■ VV rates are in excellent agreement with experiment. 
Direct CP asymmetries 

1) . In the heavy quark limit, the predicted CP asymmetries for the penguin-dominated modes K^n^, 

K*^K^, K p^, K^p^, and tree-dominated modes tt+tt", p^K^ (with A^p defined in Eq. (15.191 )) and 
p^Ti^ are wrong in signs when confronted with experiment. Their signs can be flipped into the right 
direction by the power corrections from penguin annihilation. 

2) . On the contrary, the decays K^-nP, K r\, K*^r\, 7r°;r° and n^r\ get wrong signs for their direct 

CP violation when penguin annihilation is turned on. These CP puzzles can be resolved by having 
soft corrections to the color-suppressed tree coefficient so that a2 is large and complex. 

3) . The smallness of the CP asymmetry in B n^n^ is not affected by the soft corrections under con- 

sideration. This is different from the topological quark diagram approach where the color-suppressed 
tree topology is also large and complex, but Acp{7i^7t'^) is predicted to be of order a few percent. 

4) . If the color-suppressed tree and electroweak penguin amplitudes are negligible compared to 

QCD penguins, CP asymmetry differences of K^k'^ and K^K^, K^K^ and K^^K^, K*^7i^ and 
K*-K+, K*^K^ and ^*°;r" will be expected to be small. Defining AAf.(,)^ = Acp{K^*'>^ 7i^) - 
Acp{K^*)-%+) and ^A'^^„^ =Acp{K^*^''n'') -Acp{K^*^H-), we found AAk^ = (12.3l^-0)%, aA^^ = 
(-11.0+^;^)%, AAK*n = {U.ltji)% and AA^.^ = (-11.1+^;^)%, while they are very small (less 
than 2%) in the absence of power corrections to the topological amplitude c'. Experimentally, it will 
be important to measure the last three CP asymmetry differences. 

5) . For both 5° K^n^ and B° K*^n^ decays, their CP asymmetries are predicted to be of order 

—0.10 (less than 1%) in the presence (absence) of power corrections to a2. The relation AA^jj. « 
—AAkii and the smallness of Acp{KPn^) give a model-independent statement that Acp{K^'K^) is 
roughly of order —0.15. Hence, an observation of Acp{SPn^) at the level of —(0.10 ~ 0.15) will 
give a strong support for the presence of soft corrections to c' . It is also in agreement with the value 
inferred from the CP-asymmetry sum rule, or SU(3) relation or the diagrammatical approach. For 
gO ^ ^OpO^ obtained Acp(^V) = 0-087+°;°^^. 

6) . Power conections to the color-suppressed tree amplitude is needed to improve the prediction for 

Acp{K*^r]). The current measurement Acp{K*^r]) = 0.19 ± 0.05 is in better agreement with QCDF 
than pQCD and SCET. 

7) . There are 6 modes in which direct CP asymmetries have been measured with significance above 3a: 

K-7i+,7i+K-,K-T],K*°T],K-p^ and p^Tl"^. There are also 7 channels with significance between 
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3.0(7 and 1.8a for CP violation: p^K^ ,K*^K^ ,K^K^ ,K^ri ,(oK^ ,71^71^ and p^K^ . We have shown 
in this work that the QCDF predictions of Acp for aforementioned 13 decays are in agreement with 
experiment except the decay 5° (oK^. The QCDF prediction Acp{(oK^) = -OMI^qqH is not 
consistent with the experimental average, 0.32 ±0.17. However, we notice that BaBar and Belle 
measurements of Acp{coK^^) are of opposite sign. 

Mixing-induced CP asymmetries 

a) . The decay modes rj'Ks and ^Kg appear theoretically very clean in QCDF; for these modes the central 

value of ASf as well as the uncertainties are rather small. 

b) . The QCDF approach predicts ASj^oj^^ « 0.12, ASaKg ~ 0.17, and ASpOj^^ —0.17. Soft corrections 

to 02 have significant effects on these three observables, especially the last one. 

c) . For tree-dominated modes, the predicted S^+j^- ^ —0.69 agrees well with experiment, while Spo^^o ~ 

—0.24 disagrees with the data in sign. 

Puzzles to be resolved 

i) . Both QCDF and pQCD can manage to lead to a correct sign for Acp{K^'f]), but the predicted magni- 

tude still falls short of the measurement —0.37 ±0.09. The same is also true for Acp(7r^7r^). 

ii) . The QCDF prediction for the branching fraction of B ^ K*r]' , of order 1.5 x 10^^, is smaller com- 

pared to pQCD and SCET. Moreover, although the QCDF results are smaller than the BaBar measure- 
ments, they are consistent with Belle's upper limits. It will be crucial to measure them to discriminate 
between various predictions. 

iii) . CP asymmetry of B° (oK^ is estimated to be of order -0.047. The current data 0.52+o 2o ± 0.03 by 

BaBar and -0.09 ± 0.29 ± 0.06 by Belle seem to favor a positive Acp{(oK^). This should be clarified 
by more accurate measurements. 

iv) . CP violation of p'^Tr" is predicted to be of order 0.11 by QCDF and negative by pQCD and 

SCET. The current data are 0.10 ±0.40 ±0.53 by BaBai" and -0.49 ± 0.36 ± 0.28 by Belle. This 
issue needs to be resolved. 

In this work we have collected all the pQCD and SCET predictions whenever available and made a 
detailed comparison with the QCDF results. In general, QCDF predictions for the branching fractions and 
direct CP asymmetries of B ^ PP, VP, VV decays are in good agreement with experiment except for a few 
remaining puzzles mentioned above. For the pQCD approach, predictions on the penguin-dominated VV 
modes and tree-dominated VP channels should be updated. Since the sign of Acp{K^ri) gets modified by 
the NLO effects, it appears that all pQCD calculations should be carried out systematically to the complete 
NLO (not just partial NLO) in order to have reliable estimates of CP violation. 

As for the approach of SCET, its phenomenological analysis so far is not quite successful in several 
places. For example, the predicted branching fraction 3S{B^ P V') ~ 0-4 x 10^^ is far below the ex- 
perimental value of ~ 9 X 10^^. The most serious ones are the CP asymmetries for K^7i^\ 7r"7r'', n^rj and 
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K*^ri . The predicted signs of them disagree with the data (so the AAxn puzzle is not resolved). Also the pre- 
dicted CP violation for K^n^ and K*^7i^ is of opposite sign to QCDF and pQCD. As stressed before, all the 
B-CP puzzles occun^ed in QCDF will also manifest in SCET because the long-distance charming penguins 
in the latter mimic the penguin annihilation effects in the former. This means that one needs other large and 
complex power corrections to resolve the CP puzzles induced by charming penguins. For example, in the 
cun^ent phenomenological analysis of SCET, the ratio of C'^ /T^'^ is small and real to the leading order. This 
constraint should be released somehow. 
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